KEK-TH-344 

KEK preprint 92-86 

CHIBA-EP-65 

DPNU-92-23 

October 1992 



Renormalization in the Gauged Nambu-Jona-Lasinio 

Model 

Kei-ichi Kondo* 
Department of Physics, Chiba University 
Chiba 263, Japan 

MASAHARU TANABASfflt* 
National Laboratory for High Energy Physics (KEK) 
Tsukuba, Ibaraki 305, Japan 

Koichi Yamawaki § 
Department of Physics, Nagoya University 
Nagoya 464-01, Japan 



Submitted to Prog. Theor. Phys. 



*E-mail address: kondo@tansei.cc.u-tokyo.ac.jp 

t Fellow of the Japan Society for the Promotion of Science. 

^E-mail address: tanabasMtheory.kek.jp 

§ E-mail address: b42060a@nucc.cc.nagoya-u.ac.jp 



1 



Abstract 



Based on the Cornwall-Jackiw-Tomboulis effective potentiai, we extensively 
study nonperturbative renormalization of the gauged Nambu- Jona-Lasinio model 
in the ladder approximation with standing gauge coupling. Although the pure 
Nambu-Jona-Lasinio model is not renormalizable, presence of the gauge interac- 
tion makes it possible that the theory is renormalized as an interacting continuum 
theory at the critical line in the ladder approximation. Extra higher dimensional 
operators ("counter terms") are not needed for the theory to be renormalized. 
By virtue of the effective potential approach, the renormalization ("symmetric 
renormalization") is performed in a phase-independent manner both for the sym- 
metric and the spontaneously broken phases of the chiral symmetry. We explicitly 
obtain (3 function having a nontrivial ultraviolet fixed line for the renormalized 
coupling as well as the bare one. In both phases the anomalous dimension is 
very large (> 1) without discontinuity across the fixed line. Operator product 
expansion is explicitly constructed, which is consistent with the large anoma- 
lous dimension owing to the appearance of the nontrivial extra power behavior 
in the Wilson coefficient for the unit operator. The symmetric renormalization 
breaks down at the critical gauge coupling, which is cured by the generalized 
renormalization scheme ("M-dependent renormalization"). Also emphasized is 
the formal resemblance to the four-fermion theory in less than four dimensions 
which is renormalizable in 1/N expansion. 
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1 Introduction 



This is an expanded version of our previous paper [JT| on the renormalization of the 
gauged Nambu-Jona-Lasinio (NJL) models (gauge theories plus NJL-type p[ four- 
fermion interactions). 

The gauged NJL models have recently become very popular in the context of mod- 
ern versions of the dynamical electroweak symmetry breaking such as the walking 
technicolor [Q, technicolor models with strong coupling ETC|5], ||, top quark conden- 
sate model (top mode standard model) [0], etc., and also in the context of a possible 
existence of nontrivial (interacting) QED[|| [J. Among others the most important 
feature of the gauged NJL model is a very large anomalous dimension [p~0| 



l< 7m <2, (1.1) 
which corresponds to a very slowly damping behavior of the fermion dynamical mass, 

£(-p 2 ) ~ (-pT 1+W2 - (1-2) 

The simplest version of the gauged NJL model, quenched QED plus chiral invari- 
ant four-fermion interaction {G / 2)[{ipip) 2 + (ipi^ip) 2 }, was first studied by Bardeen, 



Leung and Love|TT]] in the ladder Schwinger- Dyson (SD) equation. A full set of spon- 
taneous chiral symmetry breaking (S^SB) solutions of the ladder SD equation and the 
critical line were discovered by Kondo, Mino and Yamawaki[|l2]| and independently by 
Appelquist, Soldate, Takeuchi and Wijewardhana[|l3l . The critical line reads(Fig. P; 

= g* (0<a<a c = ^), (1.3) 

a = a c ='- (g<^), (1.4) 

where a = e 2 /4ir and g = GA 2 /4n 2 , with A being the ultraviolet cutoff. This is the 
line separating the S^SB phase (g > g*) and the unbroken (symmetric) phase (g < g*). 
The critical line Eq. (|1.4j) is actually the nontrivial ultraviolet (UV) fixed line [12, 




|14| , |I5| , with the line of a = constant being identified as the renormalization-group 
(RG) flow. This identification is in accord with the usual expectation that the gauge 
coupling a may not be renormalized in the absence of the vacuum polarization in the 
ladder approximation. Actually, this model, having non-running gauge coupling, may 
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be regarded fl~6|| as the "standing" (non-running) limit of the "walking" (slowly running) 
gauge theories plus four-fermion interaction in the "improved" ladder approximation 
(ladder SD equation with the gauge coupling simply replaced by the one-loop running 
one)|g, [17], g]§ @§.|] Once the RG flow is so identified, the scaling relation^ [13 



M d g-g*\ iV^T^c 

(g>g) (1.5) 



A \g-g*J 
implies an explicit form of the f3 function for |21 



P g {g, a) = A^\ a ,M d = -2(5 - <?*)(<? - ?) (g > g*), (1-6) 



with M d = S(0) and g* = | (l — yl — a/a c j . Eq.(|i~6|) indeed has a nontrivial UV 
fixed line at g = g*. 

The very large anomalous dimension was in fact found at the UV fixed line by 



Miransky and Yamawaki RIO 



7m=i + */i-- (g = g*), (i-7) 

V a. 



which corresponds to the slowly damping S^SB solutions obtained by Refs.[|12], 13 



S(-p ) ~ -P v ■ (1.8) 



It was further suggested by Miransky and Yamawaki [pH|1 that such a large anoma- 
lous dimension 7 m > 1 would imply the (nonperturbative) renormalizability of the 
four-fermion interaction, since the four-fermion operators would then become rele- 
vant/marginal, d(jj^\2 = 2d^ = 2(3 — 7 m ) < 4, in the ladder approximation. More- 
over, due to the presence of gauge coupling a ^ (j m < 2) the theory might have a 
nontrivial (interacting) continuum limit A — > oo in contrast to the pure NJL model 
with a = ( 7m = 2)© 0,0. 

However, the above nonperturbative renormalization procedure, originally proposed 
by Miransky || in the ladder QED (without four-fermion interaction), has so far been 
made only for the bare couplings and for the S\SB phase. We wish to find the RG 
property in terms of the renormalized couplings and in the symmetric phase as well. 



1 For the effects of the vacuum polarization in QED plus four-fermion interaction see Ref . | 
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Indeed, running of the bare couplings based on such a renormalization in the restricted 
coupling space may not correspond to the conventional (3 function of the renormalized 
couplings in the continuum theory, unless the RG flow is correctly identified. 

Furthermore, the above renormalization procedure was crucially based on the non- 
trivial (S^SB) solution of the SD gap equation for the fermion mass function, which 
is no longer possible in the symmetric phase where the gap equation has only a trivial 
solution. This would yield an identically vanishing (3 function for the bare couplings 
in the symmetric phase (non-running for g as well as a). Accordingly, the anomalous 
dimension in the symmetric phase of the gauged NJL model was considered to be small 
(7 m = 1 — 4/1 — — < 1)|P-1]| even in the vicinity of the UV fixed line, which is contrasted 



with that of the S^SB phase (7™ = 1 + .^/l — ^ > 1) |JTp[] , thus implying a paradoxical 
discontinuity of the anomalous dimension across the UV fixed line. 



However, it was pointed out by Kikukawa and Yamawaki|23| that such a disconti- 
nuity would be an artifact of non-running treatment of the four-fermion coupling g in 
the symmetric phase.Q A possible resolution was in fact demonstrated |23f in the four- 
fermion theory in D (2 < D < 4) dimensions (to be generically denoted by NJLd<4 
hereafter) where nonperturbative renormalization can be explicitly done in the 1/N 
expansion [ p5|1 . Through the renormalization of the fermion four-point function (aux- 
iliary field propagator) as well as the two-point function (fermion propagator), one 
obtains running of the coupling in the symmetric phase as well as the S^SB phase. 
The /3 function does have a nontrivial UV fixed point not only for the bare coupling 
but for the renormalized coupling. This in fact gives rise to a large anomalous dimen- 
sion ( / y m = D — 2) near the UV fixed point, thus filling in the would-be discontinuity of 
the anomalous dimension across the UV fixed point in NJLb<4[^|. The large anoma- 
lous dimension does exist not only for the bare coupling but also for the renormalized 
coupling in the continuum theory. Such a large anomalous dimension was in fact ex- 
plicitly shown |23|] to be consistent with the operator product expansion (OPE). Most 



remarkably, the Wilson coefficient for the unit operator does have an extra nontrivial 
power behavior (other than the anomalous dimension) due to nonperturbative effects. 

In the previous paper[I|] we showed that thanks to the presence of gauge interactions 
(a 0), the gauged NJL model can be renormalized in the ladder approximation in a 
quite similar fashion to NJLo<4- The crucial point was that the renormalization was 
done through the effective potential in the symmetric as well as the S\SB phase in a 



2 The nontrivial scaling behavior near the critical line in the symmetric phase was also suggested 
through the gap equation for the pure NJL model in Ref. p4| . 
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phase-independent manner, in contrast to the earlier works based on the SD gap equa- 
tion. To demonstrate such an advantage of the effective potential approach, we first 
reformulated the renormalization procedure of Ref.pBfl for NJLd<4 through the effective 
potential (Similar reformulation for NJLd<4 was also made by Ref.[Q). Then, for the 
gauged NJL model in four dimensions we considered the Cornwall- Jackiw-Tomboulis 
(CJT) effective potential[^7f and rewrote it only in terms of the local auxiliary fields 
a la Bardeen and Love[p28[. This effective potential, an analogue of the effective po- 
tential for NJL£><4, was then renormalized in a very similar manner to NJLo<4 and 
was explicitly written in terms of the renormalized parameters of the continuum limit 
theory (A — > oo). Remarkably enough, as in NJLd <4 the auxiliary field propagator |29 
was shown to be simultaneously renormalized through the above renormalization of the 
effective potential. We explicitly computed the (3 function for g for the renormalized as 
well as the bare coupling, which in either case has a nontrivial UV fixed point for each 
a (fixed line in (ct,g) plane). In either case we obtained a large anomalous dimension 
both in the SxSB and the symmetric phases without paradoxical discontinuity across 
the UV fixed line, in accord with Ref . ||23|| . As in NJLd<4 the OPE was explicitly given 
in a consistent manner with such a large anomalous dimension in both phases. As in 
NJL £) <4 the Wilson coefficient for the unit operator acquires an extra nontrivial power 
behavior other than the anomalous dimension. 

In this paper we present detailed description of the results of Ref. jT]] on the renor- 
malization of the simplest gauged NJL model, gauge theories with standing gauge 
coupling plus four-fermion interaction, in the ladder approximation. As a basis of our 
analysis we consider the CJT effective potential written in terms of the auxiliary fields. 
Here we give a more general form than the simplest one (the BL form|28|]) discussed 
in Ref. |l|]. By use of this effective potential, it is shown in the S%SB phase that all 
the amputated multi-fermion Green functions at zero momentum are finite at the crit- 
ical line (including the end point a = a c ) in the continuum limit A — > oo. We then 
give an explicit procedure to renormalize this effective potential in that limit. This 
renormalization is possible owing to the presence of gauge interactions fa / Oj. The 
(3 function and the large anomalous dimension are obtained through this renormaliza- 
tion both in the symmetric and the SxSB phases. The OPE is explicitly constructed, 
which is consistent with the large anomalous dimension in both phases. Correspond- 
ing to the generalized form of the effective potential, we consider a generalization 
( "M- dependent renormalization") of the simplest renormalization scheme ("symmet- 
ric renormalization" ) [jlj made on the symmetric vacuum. Wilson coefficients and RG 
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functions, etc. are calculated in the M-dependent renormalization as well as in the sym- 
metric renormalization. In particular, whereas the symmetric renormalization breaks 
down at the end point a = a c of the critical line, the M-dependent renormalization 
still remains valid there. 

The paper is organized as follows. In the next section we derive the CJT effective 
potential and its variants of the gauged NJL model (in the equivalent Yukawa form 
rewritten in terms of the local auxiliary fields) in the ladder approximation. In section 3 
all the amputated multi-fermion Green functions at zero momentum are explicitly cal- 
culated from the effective potential in the S^SB phase and shown to be finite even at 
the end point a = a c of the critical line in the continuum limit A — > oo. Section 4 
critically reviews calculation of the auxiliary field propagator made by Appelquist et 
al.[^]. Then in section 5 we present an explicit procedure of renormalization (sym- 
metric renormalization done on the symmetric vacuum) which is made through the 



effective potential (BL effective potential p8[) in an analogous manner to the renor- 



malization of NJL£> <4 |23|]. In section 6 explicit construction of OPE is given, which 
is shown to be consistent with the large anomalous dimension in both the symmet- 
ric and the S^SB phases in a quite nontrivial manner: The Wilson coefficient for the 
unit operator possesses an extra power behavior other than the anomalous dimen- 
sion. Section 7 is devoted to the M-dependent renormalization, a generalization of the 
symmetric renormalization, which remains valid at the end point a = a c where the 
symmetric renormalization breaks down. Section 8 is the conclusion and discussion: 
We comment on the "renormalizabilty" of NJL£><4 and the gauged NJL model in the 
language of the usual RG equation of the equivalent Yukawa model. In Appendix A 
the effective potential and auxiliary field propagators in NJL D<4 are given. Appendix 
B presents OPE in NJLo<4. In Appendices C and D we present RG study of the bare 
parameters d la Miransky through the SD equation and through the effective potential, 
respectively. 

2 Effective Potentials of Gauged NJL model 

Let us start with the lagrangian of the SU(N) gauge theory plus NJL-type four-fermion 
interaction: 

£ = - em - moijij + ^ \{^f + (^75^) 2 1 - ^tr(F^F^), (2.1) 



2 
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where mo is the bare fermion mass, e the gauge coupling constant, and G the four- 
fermion coupling. By using auxiliary fields a, it, Eq.(2.1) is cast into an equivalent 
lagrangian 



C = ${0 - em ~ + ™7 5 )^ - V {d) (a, tt) - -tr(F^F^), (2.2) 
where the classical part of the potential of auxiliary fields a and ir is given by 



V( c i)(o-,7t) 



1 

G 



N 



a 2 + 7T 2 I — m cr 



(2.3) 



It is more convenient to study Eq^^) than Eq. (|2.1| ) for the discussions of renormal- 
ization. 

The point is that as is demonstrated in NJLd<4[[1], the renormalization is studied 
most transparently through the effective action (potential) written only in terms of 
the auxiliary fields. Such an effective potential in the gauged NJL model is derived 
from Eq. Q2.2p by integrating out the degrees of freedom of fermion (-0) and gauge boson 
(An) . An effective potential of this kind was first derived by Bardeen and Love (BL) p8| 
through discussion of the vacuum condensate of fermion composite operator. Here 
we take an alternative approach[||] using the effective action of Cornwall, Jackiw and 
Tomboulis (CJT)|27]. [3(J and derive systematically several variants of the CJT effective 
potential including the BL effective potential as a special case.f] We thus clarify the 
relation among various kinds of effective potentials including the BL potential. 

In the CJT formalism IE7H we introduce a bilocal external source in a similar manner 



to the usual external source terms: 



Aource = -ip{x) J(x, y)ip(y) , (2.4) 

in the generating functional: 

W[J;a,n] = -2 In J [#][#] [gauge] exp [i J d 4 x(£ + £ SOUICC )]. (2.5) 

Corresponding to the bilocal external source term J, the fermion propagator S becomes 
3 It was also noted in Ref. |55) that the BL potential can be derived from the CJT effective potential. 
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a variational variable in the CJT effective action through the Legendre transformation, 



T[S; a, tt] = W[J; a, tt] - Tr( J • 5) 

= -«Tr (hnS' 1 + S^S) + k 2PI [S] - / d*xV {cl) (a, tt), (2.6) 

where k 2PI [5'] is the sum of all the two-particle (fermion) irreducible diagrams written 
in terms of S, and So is the function of a, 

iSg 1 = — a — 7TZ7 5 . (2.7) 

The stationary condition for S gives the SD equation of the fermion propagator 

5T , , 5k 2PI 

= J = z— = -S- 1 + S^ + i— — . (2.8) 

<>s ss 



Since the SD equation Eq.(2.£) can be regarded as the condition for the bilocal external 



source to vanish, Eq.(|2.6j) reads 

T[S sol ;a,n] = W[J = 0;a,n], (2.9) 



with S so \ being the solution of Eq.( |2.8| ). This in fact yields the desired effective action. 
Note that S'soi depends on the auxiliary fields o and tt. 
The Yukawa-type vertex Ts is given by: 

Ts(x,y;z) = -ij^jS~l(x,y), (2.10) 

which satisfies the following SD equation: 



T s (x,y,z) = 5^(x-y)5^(y-z) 
i5 2 K 2m [S] 



+ 



5S(y,x)5S(x',y' 



S so i(x', x")T s (x", y"; z)S so i(y", y'), (2.11) 



where integration over the repeated indices x',x",y',y" is understood. 

Now, noting the translational invariance, S(x,y) = S(x — y), cr(x) = a and tt(x) = 
tt, we define the CJT effective potential V; 

V[S; a, tt) = —T[S, a = const, tt = const] /fi, (2-12) 
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with Q being the space-time volume. Then Eq. (|2.12 ) may be rewritten as 



V[S; a, tt) = V (d) (a, vr) + V {qn) [S; a, vr), (2.13) 
where the quantum part V( qu ) [S; a, tt) is given by 

- ttV (qu) [S; a, tt) = -^^(LnS- 1 + S^S) + n 2Pl [S]. (2.14) 

2.1 Ladder approximation 

For actual calculation we need to make an approximation for K^fS] which contains 
infinite number of diagrams. Here we consider the simplest choice, namely, the lowest 
diagram (two- loop diagram) depicted in Fig. |2] 4 . This corresponds to the ladder 
approximation with the fixed gauge coupling. Although the ladder approximation is 
not a systematic expansion, it picks up at least the leading log behavior of fermion 
mass function which is actually important for our present purpose to renormalize the 
gauged NJL model. Since the running effects of the gauge coupling is left out of 
account in this approximation, it is certainly not a good approximation for the QCD- 
like gauged NJL model with a normal running gauge coupling. However, it can be 



regarded |16| as the standing limit of the walking gauge theories (gauge theories with a 
slowly running coupling) plus four-fermion interactions. For N = 1 this approximation 
also corresponds to the quenched QED plus NJL-type four-fermion interaction. 

In the ladder approximation we may parameterize the fermion propagator S (in 
Landau gauge) by 

zS-'ip) =i>- S(-p 2 ) - *7 5 £ 5 (-p 2 )- (2.15) 
Then we evaluate each term of the quantum part of the CJT effective potential 



TrLnS" 1 /-A d 4 p 

-trlnS 1 (p) 



iQN J (2tt 



A 2 



£V £ ) + £^ 2 ' 



d P y E Hl + y " E ' - b ^ EJ ), (2.16) 



lir 2 Jo pi 



4 Our effective potential is denned so as to keep the auxiliary fields not path-integrated out and 
hence does not include the "two-loop" graph with the a, tt line sitting on the diameter of the fermion 
line circle. Or, even if we included such a diagram, it would vanish identically anyway, because u, tt 
are not propagating at this stage. 
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d 4 p 



■tr [{$ 



a 



ij 5 7T)S(p)] 



Aix z Jo 



A 2 



i 2 2 

dp E p l 



<jZ( P 2 B )+TrZ 5 {p 



p| + E 2 (p|) + ^( F 



£ 2 (p 2 )+£ 2 (p : 



pi + H\pl) + E 2 (p 



,(2-17) 



a d 4 p d 4 k 



(2%) 4 i (2vr) 4 z 



tr[-fS(pWS(k)]Dr,(p-k), 



8tt 2 



A 2 



dp\p\ 



Qjf\i r-, rv i-, 



(p 2 +£ 2 + £ 2 )(A; 2 +£ 2 + £ 



K(p 2 ,^), (2.18) 



where the Euclidean momentum integral (p 2 E = — p 2 ) is regularized by the ultraviolet 
cutoff A, Cf is the quadratic Casimir of the fermion representation, and the gauge 
boson propagator D^ v in Landau gauge takes the form D^q) = (—i/q 2 )(g^—q^qu/q 2 ), 
and 

2 ,2\ SC F a/A7r 



K(p 2 E , K 



(2.19) 



max(p 2 , k 2 E ) 

We have defined our effective potential by subtracting a variable-independent diver- 
gence from Eq. (|2.14|) at the origin E = £ 5 = 0, o = n = 0: V( qu )[£ = 0, E 5 = 0; cr = 
0,tt = 0) = 0. 

Plugging Eqs.( p.l6| - p.l8D into Eq. ( |2.14| ), we obtain 



An 2 



N 



■V[£,£ 5 ;0-,7r) 



A 



2 r 



(a 2 + 7T 2 ) — m CT 



A 2 



+ 



+ - 

2 



dp 2 p 2 



In 1 + 



£ 2 + E^ 



S 2 + E 2 



(jE + 7rEc: 



pI 



+ 



p 2 + E 2 + E 2 p 2 + E 2 + E 2 



1 /-A 2 



rfp|p 2 



A 2 



EG^E^ + E^ 2 )^ 2 



dktkt-—- - — - —K(p 2 E , fc|), 

(p 2 +E 2 + E 2 )(A ; 2 +E 2 + E 2 ) 



(2.20) 



with 



GA 2 
4^ 



(2.21) 
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where we included the classical part Eg. (|2.3|). It should be noted that our cutoff 
regularization does not violate the chiral symmetry: 



cos#, sin 8 
— sin 9, cos 9 













) 


\ " ) 


• 







cos#, sin 6 \ I E 
— sin#, cos# / \ E 5 



(2.22) 



2.2 V[E,£ 5 ] 



Starting with the CJT effective potential Eq. (|2.20| ), we now investigate its variants in 
what follows. We first derive a form of the CJT potential written only in terms of 
dynamical mass E and E 5 , by eliminating a and tt through their stationary conditions. 
The stationary condition of the effective potential for auxiliary fields gives 



4tt 2 d 
N da 



V[E, S 5 ;cr,7r) 



A 2 p^E 
" ^pI + S 2 + Eg 



4vr 2 d 

= ^[E,E 5 ;a,7r) 



N dir 

t-A 2 



dpi 

" Pi + s a + s 5 



A 2 

+ -(* 
9 



A 2 

+ — 7T. 



m J 



(2.23a) 



(2.23b) 



Plugging the solution of Eq.( p.23| ) back into V[E, E 5 ; a, 7r), we obtain the effective 
potential written only in terms of the mass function of the fermion: 



4tt 2 
N 



■F[E,E E 



o 



dplpl 



1 E 2 
-ln(l + — 
2 



E 2 



E 2 

^5 



P 2 , 



+ 



m E 



A 2 



+ - / ^P 2 P 2 

2 Jo 



p| + E^ + E^ p 2 +E^ + E^ 



A2 2 , 2 S(p 2 )S(A; 2 ) + S 5 (p 2 )S 5 (A; 2 ) 

LLrC r-,rC ^ 

; - - ■ 2 v2 , ;,2 v2 



(p 2 +E 2 + E 2 )(fc 2 +E 2 + E 2 



A 2 



.(2.24) 



This type of CJT potential was obtained by Nonoyama, Suzuki and Yamawaki||l4 



directly from the original lagrangian Eq. ( |2.1| ) with the lowest k 2P1 being given by Fig. |3|. 
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Actually, the stationary condition of this effective potential Eq. ( |2.24| ) leads to the usual 
ladder SD gap equation for the fermion mass function Eq. (|C.6|) in Appendix C. 



2.3 ^[E sol ,S s 5 ol ;(j 3 7r) 

We now derive explicit form of the effective potential written solely in terms of the 
auxiliary fields through several intermediate steps described in this and the next two 
subsections. 

As such an intermediate step we first obtain a variant of the CJT effective potential 
l/[E so1 , E| o1 ; a, 7i") by plugging the solution E so1 and E5 01 of the stationary conditions, 
5V/5T, = 0, SV/SHs = 0, back into V[E, E 5 , a, ir). The stationary conditions (still not 
"gap equations") are equivalent to the ladder SD equations 



s 5 (p 



dkl klnkl 

B fc 2 +E 2 + E 2 
'0 ^ 2 + E 2 + E 2 



2 U2\ 
E)1 



K(p 2 E ,k 
K(P 2 E , k 2 E 



(2.25a) 
(2.25b) 



which are reduced to Eq.( |C.5| ), the stationary condition of Eq.( [2.24|) , when a and 7r 
are eliminated by use of Eq.( 2.23| ). 

By using the chiral symmetry, we can always rotate E,E 5 ,<7 and it so as to set 
E5 = 7T = 0, where E and E5 are understood to be the solution of Eqs. (|2.25a| - |2.25b| ). 
Thus, it is sufficient to study the case of E 5 = and n = 0. The SD equation Eq. fl2.25a]) 
now reads 

,2 



= &+ I dk E 

Jo kl + E 2 



(2.26) 



For actual evaluation of y[E sol ,E 5 = 0; a, n = 0) it is useful to note|36|, [b|] that 
V[E, 0, a, 0) obeys a simple scaling relation: 



V (qu) [S K , 0; a, 0; A) = « 4 V (qu) [E, 0; <r/«, 0; A/ 



(2.27) 



with E K (p|) = kE(p|/k ), where we made explicit the A-dependence of Vr qu ). Taking 
k derivative of Eq. ( |2.27| ) at k — 1, we find 



V (qu) [E sol ,0;a,0;A) = ^ 



2A 2 



d 
dA 2 



a- 



d 

da 



^ (qu) [E sol ,0;a,0;A), 



(2.28) 
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where we have used <5V( qu )/<5£| 2=S soi = 0. By using quantum part of Eq. ( |2.20| ) and 



Eq.( 2.28 ), we obtain 
4tt 2 



\/ (qu) [S sol ,0;a,0;A) 



N 

A 4 I 1 



, In I 1 + — 
2 I 2 I A 2 



^A 



crE/ 



A- 4 EH + A 2 + E^ 



3C F A 2 E, 



-a- 



8tt A 2 + E 2 , -/o 



A 2 



dk 



E kl + (E so1 ) 2 



+ 



fc 2 S sol (A; 2 ) cr M 2 2 p 2 £ S sol (p 2 ' 



dp' 

4^o % 2 + (E s 



where 



E A = S sol (p 2 = A 2 ). 
Putting p| = A 2 in the SD equation Eq. (|2.26| ), we obtain 



(2.29) 



(2.30) 



E A - cr 



3CW f A2 dk 2_^E^_ 
4tt A 2 -A) B fc 2 + (E so1 ) 2 ' 



(2.31) 



Plugging Eq.( p.31|) into Eq. (|2.29|) , we obtain another expression for the CJT effective 
potential 



An' 



^[E sol ,0;a,0;A) 



N 



A 



2 r 



-cr — m cr 



y2 



+ ^-ln 11 + —) ! 

I I A / 3Cf ot 



4tt A 2 (E a -ct)ct 



(2.32) 



where the classical part Eq. (|2.3|) was included. Note that Ea is a function of cr as 
determined by the SD equation Eq. Q2.26Q and hence Eq. ( |2.32| ) can in principle be 
written only in terms of cr. 



2.4 Fermion mass function 

In order to solve Ea as an explicit function of cr, we now discuss the solution E so1 of 
the SD equation Eq. ( |2.26| ) fli~2| , |13|]. Eq.( [2.26| ) is equivalent to the differential equation 



pI 



d 

dpi, 



+ 2 



+ 



3CV 



a 



dpi 4vr p 2 +E 2 ( P 2 ; 



E( F 



(2.33) 



14 



with infrared (IR) boundary condition (BC): 



d 

Pa- dp 2 E 



(2.34) 



and ultraviolet (UV) BC: 



pI 



d 

dp 2 E 



a. 



(2.35) 



In high energy region p 2 E ^> S 2 (p|), the differential equation Eq. ( |2.33|) can be safely 
approximated by the linearized equation 



Pi 



d 

~l 



d 3C F a 
+ 2— + 



dp E I dp 



^ Pi 



npi) + o( 



s 3 (pI) 



(p 



2 \2 
e) 



(2.36) 



Then the solution of Eq. (|2.36|) is written by the linear combination of two independent 
solutions: 



+d \^) +o #) <2 - 37) 

where M is an infrared scale of fermion mass function and 



CO 



a. 



(2.38) 



with 



7T 



a, = 



3CV 



The UVBC Eq. Q2.35|) determines the scale M as a function of a: 

2-w 



a ( M 

- = cA — 

A I A 



2+u) 



M 
A 



(2.39) 



(2.40) 
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where C\ and D\ are defined by 





(2.41) 



For strong gauge coupling region a > a c , u becomes pure imaginary and the fermion 
mass function becomes oscillating: 



M 



A' 



M 



sin 



/ 2 

— In -^7 + u'S' 
2 M 2 



+ 0(\ 



pl ' 

M 2 



-5/2 



(2.42) 



where 



uj = 



(2.43) 



and 



A' = 2uj'Jc 1 d l 



8' 



2iu' 



In 



(2.44) 



The coefficients c\ and d\ are complex conjugate to each other so as to guarantee that 
the fermion mass is real. The mass scale M defined in Eq. ( |2.4Uj ) becomes multivalued 
function of a for a > a c due to the oscillating behavior of the fermion mass function. 
We then take M with the largest absolute value (no- node solution), since it minimizes 
the effective potential. Note that the SD equation has a nontrivial solution even for 
cr = in this region [[31]. 

In the weak gauge coupling region < a < a c , on the other hand, the fermion mass 
function may be written as 



M 



A 



\ 



M 



sinh 



OJ Pe 

— In h u5 

2 M 



2 v -3(l-w)/2-l 



(2.45) 



where 



A = 2ujy—cidi 



5 



2(j \ d x 



(2.46) 



The coefficients C\ and d\ in Eq. (|2.37 ) are determined by the IRBC Eq.( [2.34j ). 
However, the non-linearity in the infrared region makes it difficult to calculate them 
in an analytical method. In the following, we evaluate C\ and d\ by using linearizing 
techniques of the ladder SD equation. Although the result varies slightly according to 
the choice of such a linearizing technique, we will find in sections [S], |5] and |7] that the 
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structure of the renormalization does not depend on such a detail of c\ and d\. 

There exist two familiar linearizing techniques of the ladder SD equation. One is 
to replace the £(p|) in the denominator of Eq.(OH) by M W\ M, ET|: 



(7 + 



A 2 



dk 2 K(n 2 p\ k l^( k D 
a,K B J\{p B ,K E ) ^ 2 +M 2- 



(2.47) 



Such a linearization leads to the solution: 



M 



„.l u; 1 c<j 
EY_ J 9 

l 2 2' 2 2 ' ' M 2 



oo 
n=l 



- 2 v -(l-w)/2+l-n 

M 2 ) 



+ E4 

n=l 



' p 2 \ -(l+w)/2+l-n 
M 2 J 



(2.48) 



with 



i-ir-'viuj) r ( i + | + n- i)r(-i + 1 + n - i)r(-w + i) 



dm 



. 1 Cc> .3 Co> 

r - + - r - + - 

v 2 2 ; v 2 2' 
c n (u -> -a;). 



r(i + |)r(-i + |)r ( -. + n) 



(2.49) 



In particular, we find 



Cl 



which lead to 



i r(i + cu) 

^(i + |)r(M: 



.4 



8cj cot(fw) 



di 



r(i-w) 



W r( l_« 3_u, 
v 2 2 ; v 2 2' 



\J 7T(1-^ 2 ) 

for < a < a c and 



5 = — In 

2a; 



A' 



W coth(fw') 



5' 



2iw' 



In 



1 3 cj. 

r(l + cu) r (2" 2^ ^2~ 2^ 

r(i-^)r(i + ^)r(^ + ^) 

v 2 2 ; v 2 2 ; 



r(i + ^') r (o ~^) r (o ~-o-> 



r(i-«0r(i + !£)r(§ + 2£) 



(2.50) 



(2.51) 



(2.52) 



\ 7r(l+CU' 

for a > a c . Here we took a normalization S(p| = 0) = M. It should be noted that 
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this linearizing method overestimates the order of error in Eq.( |2.37D . 

Another method is the bifurcation technique^, As a result of the bifurcation 
theory of the non-linear integral equation, the bifurcation solution from the trivial one 
satisfies the integral equation which is obtained by ignoring the S in the denominator 
and placing the infrared cutoff M in Eq. (|2.26|) : 



£(p|)=<7+ f A dk 2 E K(p 2 E ,kl)E(k 
Jm 2 



(2.53) 



The normalization of the solution is given by = M 2 ) = M. By using the bifurca- 
tion technique, we obtain 



Cl 



1 + uj 
2uj '' 



1-oj 
2u 



(2.54) 



which leads to 



10 



tanh 1 uj, (0 < a < a c ), 



(2.55a) 



A' = VTT 
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5' = — tan co', (a>a c ). 



(2.55b) 



Note that A, 5 evaluated by these linearizing methods are finite in u — > (a/ — > 0): 



A = lim A = -, 8 = lim S = ln4 - 1, 

oj— >0 71" iii— >0 



for the linearization Eq.(2.47), and 



(2.56a) 



Aq = lim A — 1 , 5o = hm 5=1, 

U— >0 a;— >0 



(2.56b) 



for the bifurcation method. 

At a = a c the fermion mass function can be obtained directly from Eq.( |2.36P plus 
IRBC Eq.(fT3§. It can also be obtained by taking uj -> (u/ -> 0) limit of Eq. flE35D 
(Eq.dggg)): 



M 



An 



M 2 



\ P 



1 Pf 

-In — + 5 

2 M 2 



(2.57) 



In this limit Eq. (|2.40| ) reads 



a = A /My 
A 2 VaJ 



, r , M 
1 + 5 - In — 



/M\ 6 

+ °<(a)> 



(2.58) 



18 



where we have used the following relation derived from Eq. fl2.46 ) and Eq. fl2.41| ): 



lim ojCx = - lim uDi = 4 £ , \im(d + D x ) = ^(5 + 1). (2.59) 



2.5 V(M) 



Let us now return to the CJT effective potential Eq. fl2.32| ). By using the explicit 
solution £ sol (j»!) of the SD equation Eq. ( |2.37|) , we find: 



S A M 
A I A 



2-u) 



(\ 2+u) 
f) Mt) 



(2.60) 



Plugging Eq. (|2.60|) and Eq . (|2.40|) into Eq. (|2.32|) , we obtain an effective potential solely 
expressed in terms of M: 



— 57T 



,V(M) 
NA 4 



1 1 

.9 9, 
1 1 



Ci 



v 4-2w / \ 4' 

Af\ 2 + cj / M 

A / 2 I A 



.9 9, 



\ 4+2w / \ 4' 

M\ 2-cu / M 

A I — + C x — 

A 2 I A 



+ 0( 



(M 



4(2-u>) 



2mo 

" a7 



\ 2— uj / \ 2+uj / \ 3(2— a;) 

C1 ,_ +c( 



where Ci and D\ are defined in Eq. (|2.41|) and 



g* = -(l + uj) 2 , ^ = -(1-0,) 
4 4 



(2.61) 



It should be stressed again that M is a function of a determined by Eq. (|2.40| ) . 



2.6 F0,tt) 

In this subsection we make more explicit the cx-dependence of the effective potential 
Eq. fl2.6ip . To this end it is convenient to rewrite Eq. ( |2.61| ) so as to leave M partly 
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unsolved: 



47r 2 V(a, 7i = 0) \m a f\ \\ a 2 



N 



A 1 



9 A' 



,9* 9 J 2A 2 



+ 



1 \ cr 



, M\ w 2 - (M 



J* 9* J 2A 2 



Ci 



M 

T 



t) 



n 2 



(2.62) 



where we have used Eq. (|2.40| ) . This effective potential is actually the basis for studying 
the renormalization in this paper, particularly in section 7. 

As a special case of our effective potential Eq. ( |2.62| ), we obtain the BL potential p8| 
derived through a different method, which is valid only for l>c<j>0(0<a< a c ). 
In this region Eq. ([2.4Up can be solved in a recursive way, 



A ) ~ dA Ci IciaJ +0 H A 



(2.63) 



with rj being 



rj = min(4. 



2(2 + u) 
2-uj 



Plugging Eq.(^63|) into Eq.(^62|) (or Eq.(gjflD), we obtain 



(2.64) 



~N~ A 1 ~ ^^ 2_+ \~9*~ 9J 2A 2 ~^ C ~Y~ \a) + VaJ j ' 

(2.65) 



where is defined by: 



-U) 



1 + u 



4/(2-w) 



c -(2 + .)/( 2 -.) rfi > Q 



(2.66) 



It is easy to recover the pseudoscalar auxiliary field 7r in the effective potential Eq.( [2.65| ) 



9 .99 / 9 9 \ 2/(2 — Cj) 

47T F(cr,7r) lm a I l\cr +7r z 4£ w 2 - u (a + tt 



AT A 4 9 A 2 \9* 9 2A 2 



a/a c 4 



A 2 
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'a 2 + 7r 2 \ v/ \ 



A 2 



(2.67) 



The expansion Eq. (|2.63| ) obviously breaks down at a = a c where (9((cr/A)' ?_1 ) term 
in Eq. fl2.63| ) gives the same order contribution as others. A remarkable feature of our 
expression Eq. fl2.62j ) for the effective potential is that it has a definite value in the limit 
u i (a t a c ): 



4vr 2 V(a, n = 0) lm a 1 A <x 2 
9 A 2 \q* 9 I 2A 2 



a 



3 r i M 

4 +5o " ln A 



iV 



A 



2A 



2 r 



In 



M 



(2.68) 



where we have used Eq.( [2.59D . Since M is written in terms of a through the UVBC 
at a = a c Eq. fl2.58Q , the above effective potential can be further expressed in terms of 
the auxiliary fiekkf] 



4vr 2 V(a, vr = 0) 1 ^o " 



1 1 



a' 



N 



A J 



9 A 2 



9 2A Z 



16- 



+ o 



where we have used the relation 



M 1, /A 2 \ ^„ , A 2 . 



(2.69) 



(2.70) 



derived from Eq.( |2.58| ). Actually, in the limit A ^> a this potential agrees with the 
effective potential derived by Bardeen and Love[^J for a = a c . 

Eq. (|2.62j ) is also applicable to the strong gauge coupling region a > a c by perform- 
ing the analytic continuation u = iu; 1 , u>' = J 'a/a c — 1. 

Finally, we comment on the relation of our derivation of the effective potential 
Eq. fl2.65|) to that by Bardeen and Love|28| which was derived based on the observation: 



d 

da 



V(qu)(0",7T = 0). 



(2.71) 



Actually, Eq. ( |2.71| ) is manifest in the C JT formalism, since V( qu ) (cr, it) is identified as 



This expression actually coincides with the one obtained directly from the solution at a = a c 
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V( qu ) [E so1 , E s 5 o1 , a, tt) which satisfies 



^-y (qu) [E so1 , o,t7,o) 



5 



da 



V( qu) [E so1 , 0,(7,0) + / 

J 



a 2 9 dE so V 



9a 5£(p' 



-V (qu) [E, )( T,0) 



S=S so1 



AT /-a 2 



£ sol G4) 

_2 /,, 9 , /v.sol\2 ; 



2„2 



where in the last line we have used the stationary condition 



(2.72) 



2.7 Gap equation 



V (qu) [S, 0,(7,0) 



5 



s=s so1 



5E(p : 



-^[E, 0,a,0) 



0. 



E=E ao1 



Now, we are interested in the solution a = <j so \ of the stationary condition 
dV(a, tt = 0)/da of the effective potential Eq. (^67|) : 



9 A 





4C 




)x- 


a/a c 


(x) 



+ 



(2.73) 



The instability of the symmetric vacuum (cr) = (7r) = for g > g* in the chiral 
symmetric limit mo = is manifest in the effective potential. Thus it is readily seen 
that 

' 5 2 (0 < a < a c ) (2.74) 



9 = 9 =-(1 + ^0 



is the critical line|12], |r3]| . Actually, we can read off the S^SB solution a = cr spont from 
Eq.( [2.73| ) at m = 0. The scaling at g ~ g* of the S^SB solution to Eq.( |2.73| ) is given 

by 



^"spoilt 

A 



a/a c 



-i (2-u)/2u 



, , , +•••, (2-75) 

where • • • stands for terms with higher power of (1/g* — 1/g). Eq.( [2.75|) combined with 
Eq. ( |2.63| ) leads to the scaling of dynamical mass of fermion M d = (M) : JT2], |T3 



Md 
A 



l/2u 



(2.76) 
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However, Eg. ( p. 67] ) is not valid in the u — > (a — > a c ) limit, in which case we need 
to return to Eg. ( 2.62 ) or Eg. ( [2.61 ), the form before the expansion Eg.( [2.63 ) is applied. 
The stationary condition of Eg.( 2.61| ) with itlq = leads to 



Md 
A 



( I 

d g* 



V 



9* 



I \ l/2w 

JL 
1 



H7 



(5-4cj)/2w 



(2.77) 



g) 



In the uj — > ( g* — > 1/4) limit we find the essential singularity-type scaling: [12|, [13 



A 



exp 



l + S - 



4-l/ 5 



(2.78) 



2.8 Yukawa- type vertex 



As shown in Eg. fl2.10|) , the Yukawa-type vertex is calculated by the a derivative of 
S^ 1 . Since we have already evaluated the fermion mass function S under the ansatz 
of constant a, it is easy to determine the Yukawa- type vertex at q = 0: 



T s ^p 2 ) = T s (p,q = 0) 



d_ 



(2.79) 



where the momentum assignment of Yukawa- type vertex Ysip, q) is depicted in Fig 
Plugging Eg.(^63|) into Eg.(p3?D, we find 



n-p 2 ) 

A 



„ x -fl-o;)/2 

y\ £ ^y 2+ .)/( 2 -^ 

1 + w \ A 2 i A +Ul U/ 



(2.80) 



which leads to the Yukawa-type vertex on the chiral symmetric vacuum (a) = 0: 

2X -(l-w)/2 



s(V 



-p 



cr=0 



1 + w A 2 



(2.81) 



In the u; —>■ (a — >• a c ), however, the expansion of the fermion mass function 
around the symmetric vacuum Eg. ( |2.80| ) loses its validity. Thus, we need to evaluate 
the Yukawa-type vertex without employing such an expansion around the symmetric 
vacuum at uj — > 0. Here we calculate the Yukawa- type vertex corresponding to the 
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expansion around non-zero a: 







n-v 2 



da 



(2.82) 



DM 



Combining the M derivative of a from Eq. ( |2.40|) , 



da 



(2-w)C 1 (^J +(2 + W )D 1 (^-J +0(( T ) ) 



and the M derivative of S(— p 2 ) from Eq.( |2.37j) 

2X -(l-u)/2 



<9£(-p 2 ) 2(2 -w) 



<9M 



we obtain 



1 + 



■Ci 



-p 

M 2 



2X -(l+w)/2 / 2 - 



1-U 



(2.83) 



-(5-3w)/2 



(2.84) 



r s (V;M) 



-(l-w)/2 



i 1 + W 

1 + K 

1 — uo 



-2+uJ 



M 2 



+ 0{ 



-v 

M 2 



1 + w A 2 



j + a ''t) +0( (t) : 



where is defined by 



K 



{2 + uj)D 1 



(2 - u)d ' 

As expected, Eq.( |2.85| ) is indeed applicable to uj — ► limit: 



(2.85) 
(2.86) 



T 5 (-p 2 ;M) = 2 



A 



1 + 25 - In 



il^ 2 
A 2 " 



(2.87) 



3 Finiteness of Amputated Green Functions 

In this section we calculate from the effective potential all the (amputated) multi- 
fermion Green functions at zero momentum of a in the S^SB phase and show that 
they are finite in the continuum limit A — > oo, once the four-fermion coupling g is fine 
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tuned so as to fix the fermion dynamical mass Md- This strongly suggests the existence 
of an explicit renormalization scheme to make the effective potential finite at the full 
critical line including the end point a = a c . 



3.1 Green functions for < a < a c 

Let us first calculate the second derivative of the effective potential Eq. (|2.62j ) 



cr 



d 

da 



V(a,n = 0), 



(3.1) 



so as to evaluate the auxiliary field propagator at zero-momentum. By using Eq. ( p.83j ) , 
we obtain 



4tt 2 yW(a) 
~N A 2 "" 



9* 9 \g* 9* , 



-(f) 



2uj 



2uj ' 



(3.2) 



where K is defined in Eq. Q2.86|) and we have neglected higher order in M/A. The 
solution of the gap equation in chiral limit, Eq.( j2.77p , reads 



1 1 

9* 9 



1 1 



(2 + lu) + (2-u)K 



Md 
A 



2uj ■ 



(3.3) 



Thus, Eq.flp]) reads 



a, 



spont / 



2 ^ ( ^f (±_± 

A J \g* g* 



N A 2 



1 + K( X 



M, N 2w 



/ /If \ 2w 
2 + uj + (2-cu)K\-^j 



(3.4) 



The cutoff dependence of ~ A 2 ^ -1 ^ cancels exactly that of the Yukawa-type vertex 
Eq. fl2.85|) and we obtain finite amputated four-point Green function in A/M^ — > oo 
limit: 



r 4 (p(i), -P(i),P(2), -P(2)) 



25 



n 

4tt 2 i=i 



on 


1 




2 


(-phY 


1 + w I M 2 d 



V s (-p\ 3)} M d ) 



(l-w)/2 / 2 

+ 2^ /Vo) 



1- w 



M d 2 



-(l+w)/2 

+ o{\ 



2 v-(5-3w)/2 



Ml 



N 



2uM?,K ( 1 1 



(3.5) 



2 + uj \g* g* 



To discuss higher-point Green functions (see Fig. |5|), we need to evaluate multi-a 
vertices. At zero-momentum of a, they are calculated from n = m + 2 (m > l)-th 
derivative of effective potential: 



V^(a)=(^-\ V(*,n = 0). 



(3.6) 



After a straightforward calculation using Eq.( |2.83| ), we find 



4^ 



(J 



;M\ -(2+m)(i-w) / i i \ 
f T J f - - - P m iy>, KiMj A) 2uj ) 



[(2 - u;)d 



2u 



2m+l 



(3.7) 



where P m (uj, z) is given by 



m—l 



(3.8) 



€=0 



with Q mt e(u) being certain polynomial in with degree m — £ — 1, e.g. 



fm = 11 



Y[((£ + l)u-2(£-l)) (m>2) 



(3.9) 



1=2 



Now, we are ready to evaluate the amputated 2n-point Green function (n > 3) at 
zero momentum. It is remarkable to see the cancellation of A dependence of Eq. (|2.85j) . 
Eq.(|3.4|) and Eq. fl3*77|) in the fermion amputated Green functions. Actually, we find the 
Green functions remain finite in continuum limit (A/ Ma — > oo): 



r 2 n(P(l), -P(1),P(2), ~P(2), ■ ■ ■ ,P(n), ~P(n)) 
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^ (n w nt ) n 



V {2) (cr spont ) 



r s (- P 2 {j) ,M d ) 



4vr2 
' N , 

n 

n 

3=1 



Qn-2fi{w) 



v -(l-w)/2 

-P(3)\ 



1 + uj \ M 2 d 



"POO 

1-w \ Mj 



-(l+w) /2 



0(| 



2 x -(5-3u)/2 



Ml 



(2uC\ 



^n-2 2^(2 - Uj) n 2 _^ 3n _ 4 ^ n _ 1 



(2 + cj) 



X* g * 



n-1 



.(3.10) 



3.2 Green functions at a = a c 

We next consider Green functions at a = a c , i.e., uj = 0. This limit should be taken 
carefully, since the Yukawa-type vertex Eq. ([2.870 and the effective potential depend on 
A not only in power of the cutoff but also in logarithm of the cutoff. As shown in the 
following, however, fmiteness of the amputated Green functions persists also in this 
limit, thanks to the cancellation of logarithmic dependence on A. 
Taking u -> limit of Eq.(^) and Eq.(|3T7|), we find 



4?r z 
~N 



16A 2 



a. 



spont j 



1 + 25, 



o 



4f 



2 + 25 - In 



'Ml 

A 2 



(3.11) 



and 



~N 



y(m+2) { 



/ M \ ~( m + 2 ) 
32M 2 " m — 
V A 



Ao 
2 



1 + 25 - In 



M 2 
A 2 



m+2 



Qm,o(0) 



2Q m ,i(0) 



1 + 2<5 n - In 



Ai 2 



+ ■ 



(3.12) 



respectively. Thus, it is clear that the logarithmic divergence of Eq. (|2.87|) cancels out 
that of Eq. (3-1 If) and Eq. (|3.12j ) in the amputated Green functions. Actually, by noting 
Qm,o(0) = (— 2) m ~ 1 (m — 1)!, we obtain finite Green functions in the continuum limit 
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A/M d -> oo: 



r 4 (P(i),-P(i),P(2),-P(2); 



--n 



-1 + 25 + In 



^1 W^f) 



(3.13) 



for the four-point Green function and 



r 2 nO(l), -P(1),P(2), -P(2), • • • ,P(t»), ~P(n)) 



4vr 2 \ n ' 1 (n-3)!A 2 

2 n M 3n-4 



iV 



n 



-1 +25 + ln 



-pco ' 



-2 1 



-P?,)/M; 



(3.14) 



for the 2n (n > 3)-point Green functions. 



4 Auxiliary Field Propagators 

We have shown in the previous section that all the multi-fermion Green functions 
are finite at zero momentum of o in the S^SB phase. This strongly suggests that the 
effective potential can be renormalized. We further wish to show finiteness of the Green 
functions also at non-zero momentum. This requires knowledge of the effective action 
not restricted to the effective "potential, which is, however, a far-reaching problem even 
in the ladder approximation. Directly relevant to such a problem are the propagators of 
the auxiliary fields. If the renormalization of the effective potential can simultaneously 
renormalize these propagators, it would be very promising for the validity of such a 
renormalization. This is actually the case in the NJL£) <4 model p3|, [1], which is 
renormalizable in 1/N expansion. 

Thus we are interested in the calculation of the propagators of the auxiliary fields. 
Difficulty in such a calculation resides in the lack of our knowledge on the Yukawa- 
type vertex function Ts{p, q) at non- vanishing momentum of the auxiliary fields q ^ 0. 
Recently, Appelquist, Terning and Wij ewar dhana [p9| 1 made an interesting resummation 
technique based on a further approximation (besides the ladder approximation) to 
evaluate the auxiliary field propagator without knowing Ts{p, q ^ 0). Most amazingly, 
as we will show in the next section, their a propagator is in fact renormalized through 
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our renormalization conditions of the effective potential. [ 



Here we briefly review the calculation of Ref.[f29j and discuss the validity and its 
possible modifications. Let us consider the propagator of auxiliary field a in the sym- 
metric vacuum cr so i = (see Fig.^J) : 



d 4 k 
(2n)H 



tr 



r s (k, q ) 



1 1 



V^(a = 0). (4.1) 



For q 2 = 0, it is given by the second derivative of the effective potential: 



iD~l{-q 2 = 0) = -V> so i = 0) 



NA / 1 



47r 2 \g* g 



(4.2) 



Now we consider the second derivative of Eq.( |4.1| ) on q^ (see Fig.|7])[[34], By 
using the self-consistent equation for T s (see Fig.|]), we can easily show that Fig.|7| can 
be expressed as 



1 g 

N dq^dq u 



-tf£(-<T) 



aTk 
(27r) 4 i 

d A k 
(2tt)H 



-tr 



tr 



r^*,,)ir s( *-,,-,) 

(ik Ts{k ' q) ) l Ts[k -"•- q) 



(H <-> u) 



(4.3) 



which is diagramatically depicted in Fig.|9|. In the same fashion, we find a symbolic 
expression of iV-th derivative: 



1 (1. 

N \dq 



d h ^ — > 



T s(k,q) TL T s (k - q, -q)- 



dq n 



dq l ~ n 



(4.4) 



It was then observed ]29] that the contribution from d n Ts(k, q)/dq n \ q= Q is negligible 
for n^O, which implies 



N \dq 



*D-J(-q 2 )\ q=0 



d 4 k 
(2tt) 4 z 



tr 



r, 



-k 2 ){-l) 



(4.5) 
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Thus, the resummation of the Taylor expansion around q = leads to a compact 
formula: 



d 4 k 
(2rr)H 



tr 



1. 

r 



r 5 (-fc 2 )^r 5 (-A^ 



+ constant. 



(4.6) 



After subtraction at q 2 = 0, the integral of Eq. ( [4.6|) yields finally: f29 



N 



A 2 



87T" JO 

N 



dk g,ht '^r 



x "e) 



8vr 
AT 



2 ?! 



8tt 



1 + 



- 2 J % 
T|(A 2 



2a c 

^ l A 2 ' 



-l+OJ 



(4.7) 



(4- 



Now, let us discuss the validity of Eq.( ^8|) . It was argued [p9[| that the contributions 
coming from 

ffr s (fc,g) 

1 9=0 



<9g n 



(4.9) 



are to be in higher order of a: 



d n T s (k,q) 



1 a 



dq n q ~° k n 4a, 



r*(-fc 2 ). 



(4.10) 



It was also shown numerically that these terms give smaller contributions in D a ^ 
than that from the zero-derivative of Tg, when the integral is regularized by an IR 
cutoff. 

However, we note that the dominant integral in Eq.( [4.6| ) comes from IR region 
where Tg(k,q = 0) does have a singularity which is actually absent in the original 
Ts(k,q 0). Thus we must be careful about the effect of higher derivatives which 
compensates such a singularity around k = 0. Here we evaluate the order of magnitude 
of such an effect, by considering possible modifications of T$ in Eq. (|4.6|) which could 
change the IR behavior of the integral of Eq. ( [4.6|) at q ^ 0. 

The simplest such modification of T$ would be the replacement in Eq. ( [4.7| ): 



r|(*S) - r 2 s (k 2 Mkl - ql) + TUql)9(ql - kl) 



(4.11) 
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Then we obtain 



iD£{fi)-iD£{0) 



Nq 2 E /3 



57T 



2 l-UJ 



(4.12) 



Another possible modification is the symmetric calculation under the exchange of k 
and k — q in Eq. (fL6|) : 



E I 



-N 



d 4 k 



(2tt)H 

Such a modification leads to 



tr 



r rf -*> ) ir,(-(*- 9 )»)yi ? 



+ constant. (4.13) 



iv r(-w) 

47r 2 r(2 + t,;) 



r(l + f) VJ 

p/3 w< 

\2 2' 



qlT> s {ql) + 0{ql). (4.14) 



Although these modifications are somewhat arbitrary, the results seem to suggest 
validity of the functional form of the a propagator: 



iD~l{-q 



2, 2 



a/ai( 



q 2 m-<f) - V"{a = 0), 



(4.15) 



with the precise form of £ w being subject to the details of the IR treatment. Whereas 
various modifications agree with each other near the pure NJL model a = 0, the 
deviation becomes significant at a = a c . Actually, £ w diverges at a = a c for Eq. ( |4.8[ ), 
while Eq. (|4.12| ) does not. In the next section we shall renormalize the generic form of 
Eq.(f4.15|) not restricted to Eq. (|4.8|), the original form of Ref. 



5 Symmetric Renormalization 

In this section we formulate the renormalization of the gauged NJL model, based on 
the effective potential Eq.( |2.62[ ) expanded around the symmetric vacuum ("symmetric 
renormalization"). This is done in a similar manner to the pure NJL model in D(2 < 
D < 4) dimensions (NJL£> <4 ) jj]]. We first reformulate [0, |26| the renormalization of 
NJL£> <4 [T23J at the leading order of 1/N through the effective potential approach. Then 
in the gauged NJL model we apply the same method through the effective potential 
written in terms of the auxiliary fields. 
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5.1 NJLd<4 model 

The lagrangian of NJLb<4 model is given by: 

G 



(5.1) 



where the fermion field ip belongs to the fundamental representation of SU(N), with 
the summation of SU(N) indices being understood. In Eq. ( |5.1[) , the fermion spinor is 
given by that in four dimensional space-time, so that the model possesses U(1)l X U (1)r 
symmetry for tjiq = besides the SU(N) symmetry. 

We first demonstrate the symmetric renormalization through effective potential. 
Introducing the auxiliary fields a, n in the same way as the gauged NJL model, we 
obtain the effective potential V(a, n) [|T 



(Att) d/2 T(D/2) 
ANA D 

1 m (T / 1 



V(a,n) 



1 \ a 2 + 7T 2 



2^ 2\ D ' 2 
O + 7T 



9 A 2 



,9 9, 



2A 2 



D 



D 



A 



0( 



a 



it 



A 2 



(5.2) 



with A being the ultraviolet cutoff in the loop integral, and g* and Cd are defined by 
g* = D/2 - 1 and ( D = B(D/2 - 1,3 - D/2), respectively. Note that g* -> and 
C^ 1 — > as D — ► 2, so that the divergence in the second and the third terms cancel 
each other to give a well-known logarithmic factor in the Gross- Neveu model |57 . For 
a detailed derivation of Eq.(|5.2|), see Appendix K|. 

Propagators of the auxiliary fields are calculated as Eq. (|A.26|) . In the symmetric 
vacuum (a) = they read 



= V"{a = 0,7i = 0) 



AN 



D 



{Atx) d I 2 T(D/2)2- D/2 
where ^ D = B(3 - D/2, D/2 - 1)/T(D - 1) and 



(V) 



2\D/2-l 



(5.3) 



V"(a = 0,n = 0) 



d 2 V(a, n = 0) 



da' 



ANA 



D-2 



1 1 



CT=0 



(Ati) d/2 T(D/2) \9 9* 



(5.4) 
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Eq.( |5.4j ) is negative for g > g*, implying appearance of a tachyon pole in the auxiliary 
field propagators, another signal of instability of the symmetric vacuum. 

Let us next consider the renormalization of NJL£j<4 model at the 1/N leading order. 
Due to absence of the divergence in the fermion propagator and the vertex function Fs, 
wave function renormalizations of the fermion and the auxiliary fields are not required 
at this stage: 

■0 = if>R, ctr = <?■ (5.5) 

Thus we concentrate our attention to the renormalization of the effective potential 
Eq.( |5.2|) . The divergence A D ~ 2 in Eq. (|5.2j ) can be absorbed into the redefinition of 
renormalized parameters g R ,m R : 

A WI_±U,W±_±) (5.6) 
\9 9') \9r 9'J 

and 



A D - 2 // D ~ 2 

m = m R , (5.7) 

9 9r 

where \i is the renormalization scale. These are precisely the same renormalization con- 
ditions as those imposed by Kikukawa and Yamawaki through the renormalization 
of the propagator of a. 

With the above definition of renormalized parameters Eqs. fl5.6| - |5.7|) we obtain a 
renormalized effective potential: 

(An) D ^ 2 T(D/2) 

-V R {a R ,Ti R ) 



AN fi D 




°\ + ^R 1 Cd ( Or + 71 'R 



D/2 



9R V 2 \9r 9r) 2/i 2 9- — D \ ^ 



(5.f 



where higher power terms 0{a i ) disappear in A — > oo (even at D = 2). This renor- 
malization breaks down at D — 4, which is signaled by the singularity (2 — D/2)~ 1 
in the last term of Eq. ( |5.8| ), corresponding to the logarithmic divergence of (a 2 + n 2 ) 2 
coupling in D = 4 NJL model. (We would need extra "counter term" such as the 
eight-fermion operators (tpip) 4 to "renormalize" the model in D = 4.) 

In our renormalization Eqs.( |5.6| - |5.7| ), the critical value of renormalized four-fermion 
coupling is left undetermined. Actually, it is a free parameter which corresponds to 
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varieties of renormalization schemes. In the D — > 2 limit, however, we need to define 



the singular part of g* R as a function of D: 

1 1 



9r D/2 - 1 



+ regular function of D, (5.9) 



so that we can reproduce the usual renormalization of the Gross-Neveu model |37 
Again, the choice of the regular part corresponds to the choice of renormalization 
scheme. In the following, we take the simplest choice g* R = D/2 — 1 = g* (different 
from that in Ref. TO- 

It is straightforward to get RG functions f3 g and 7 m from Eqs.(|5Tt3|--|5.7|): Q 



P g (g R ) = (D-2)g R \l- 9 -^y (5.10) 

7m (^) = (D-2)^, (5.11) 
9r 

where (3 g = fidgn/dfi and r ) m ra R = —fidm R /dfi. This agrees with Ref . [p3fl . It is 
evident that we also obtain the same form as Eqs.( |5.10| - p.ll| ) for /3 g (g) = Adg/dA and 



■j m (g) = — (A/m )<9m /<9A in the bare coupling g, with the UV fixed point g* R simply 
replaced by g*. As was pointed out by Kikukawa and Yamawaki it is clear from 
Eqs.( |5.10| - |5.11|) that NJLd<4 model has a nontrivial ultraviolet fixed point on which 
the RG functions are continuous when approached from both phases. 

It is also obvious that the propagators of a and it Eq. (|5.3| ) can also be renormalized 



by the same condition as Eq.( |5.6| ), as was originally done in Ref. [|23| : 



^ = 0,,, = 0) + g % a (^. (5,2) 



It should be stressed again that this renormalization breaks down at D = 4, signaled 
by the appearance of (2 — -f ) _1 singularity in the term of (q 2 E )^~ l ( — > ((?§) as D — > 4), 
which corresponds to the logarithmic divergence of the kinetic term of o and ir in D = 4 
NJL model (We would need higher dimensional operator (?/>■?/;) <9 M (ipip) as a "counter 



3 The next-to-leading order corrections to this result have been calculated [E6L p8l B9| 
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term" to "renormalize" the model in D = 4). 

The fact that we can renormalize the theory without higher dimensional operators 
("0"0) 4 an d d^ip^d^^ifj) at 1/N leading order simply reflects the following fact: (i/ji/j) 2 
is a relevant operator due to a large anomalous dimension 7 m = D — 2 at g R = g* R , 
i.e., dim(?/'^) 2 = 2(D — 1 — 7 TO ) = 2 < D, while the would-be "counter terms" (^ji/j) 4 
and d^tp^d^^ip) are irrelevant operators, dim(^-?/;) 4 = 4(D — 1 — j m ) = 4 > D, 
dim[<9 M (^V0^ M (' ? /'V ; )] = 2{D — j m ) = 4 > D. At D = 4, however, all these operators 
equally have dimension 4(= D) and become marginal operators. Hence they should be 
included in order to make the theory renormalizable, in which case the NJL model in 
its renormalizable version becomes identical to the Higgs- Yukawa system ("standard 



model") Eg, 41 



5.2 Symmetric renormalization of the gauged NJL model 

Now, we are ready to study the renormalization properties of the gauged NJL model, 
based on the simplest effective potential Eq. ([2.67| ) expanded around the symmetric 
vacuum (symmetric renormalization) |l] . The gauge coupling a does not get renormal- 
ized, in accord with the absence of vacuum polarization in the gauge boson propagator 
in this approximation. Such an approximation becomes realistic in the standing gauge 
theory as a limit of walking gauge theory [Q . Thus the renormalization is operative only 
on the four-fermion coupling g. Actually, this renormalization (0 < uo < 1) is done in a 
very similar manner to NJLd<4(2 < D < 4). Both cases break down in the pure NJL 
limit (D = 4,u = 1). However it should be noted that while the renormalization of 
NJL£> <4 is valid even at D = 2 where Eq.( |5.2|) still remains valid, the renormalization of 
the gauged NJL model in this scheme breaks down at u = where Eq.( |2.67| ) becomes 
no longer valid. Renormalization scheme valid also at uj = should be based on the 
general form of the effective potential Eq. (|2.62|) , which will be discussed in section 7. 



In our definition of (bare) auxiliary field, the (bare) Yukawa-type vertex T$ does 
not depend on the four-fermion coupling and vanishes in the A — >• oo limit. Thus, it 
should be renormalized via redefinition of auxiliary fields to obtain a finite interacting 
theory in the A — > oo limit: 





o-fl = - o", k r = - 7T, (5.13) 
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with fi being "renormalization point" .Q According to this definition of renormalized 
fields, the Yukawa-type vertex Eq.( 2.8l|) is renormalized: 



1-w 



r s (V) 



-p 



-(l-w)/2 



1 + LO \ 11 



(5.14) 



for the symmetric vacuum. It should be noted that this definition of the renormalized 
fields <jr, ti r simultaneously renormalizes the Yukawa- type vertex in S^SB vacuum 
Eq.(ET85D: 



P 2 ;M d ) 

l-u 



T s (-p 2 ;M d ) 



~P 



-(l-w)/2 



1 + U V 



1 + 



(l + o;)(2 + u;) Di / V 
(l-a;)(2-o;)C' 1 1 M j 



0(| 



2 \ 



(5.15) 



Then the effective potential Eq.( 2.67| ) is expressed in terms of the renormalized 
auxiliary fields: 

4tt 2 V R (a R , Ti R ) 



N /i 4 

l+w 

A \ 1 m a R 



9 jj 2 



1 \ ° 

9 



4£, 2-uj cr R + n 



2 2 N 2/(2-w) 



2^ 2 



a/a c 



,(5.16) 



where we have dropped out the contributions which vanish in the A — > oo limit. 

Now, the parallelism between NJL£) <4 and gauged NJL model is manifest. The 
effective potential Eq. (|2.67|) can be renormalized by the definition of renormalized 
four-fermion coupling 



(5.17) 



,9 9, 



k 9r 9r, 



7 Similar redefinition of the auxiliary fields was also made |2q| , with fj, being taken as the dynamical 
fermion mass [i — Md- 
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and renormalized current mass of the fermion: 



A 1+w — = At 1 +-_^. (5.18) 
9 9R 

Again, the ambiguity of g* R corresponds to the choice of renormalization scheme. Ac- 
cording to this renormalization, we find the renormalized effective potential 

47T 2 V R ((TR, TTr) 



N n 



1 



_2 

1 m R a R II 1 \ o% + tc r ACu, 2 - u ( o\ + nj X ~ 



• I— — 1 " " I — ' R 1 . (5.19) 

9r /i 2 \g* R 9r) 2/j 2 a/a c 4 \ /i 2 

Thus, all the multi-fermion Green functions have been renormalized at zero momentum 
of the auxiliary field. 

How about the renormalization at non-zero momentum, then? Remarkably enough, 
the auxiliary field propagators Eq.( |4.15| ) are also renormalized via the above definition 
of the renormalized parameters: 

Nf 

= V R \a R = 0, ix r = 0) + -^(-g 2 )rf \-q 2 ). (5.20) 

a/a c 

Thus the four-fermion Green function can also be renormalized at non-zero momentum 
in a similar manner to NJL£> <4 . In spite of the formal resemblance, however, one 
should note that the gauged NJL model has only been shown to be renormalized at 
the level of the ladder approximation, in sharp contrast to NJLd<4 which is shown to 
be renormalizable in the systematic 1/N expansion |25[ . 



We should also emphasize that the renormalization Eqs. ( |5.17| - |5TT^ ) is based on the 



effective potential and thus it holds both in the chiral symmetric (g R < g* R ) and the 
SxSB (g R > g* R ) phases. 

The definition of the renormalized parameters gR and tur Eqs.( |5.17| - |5.18| ) leads to 
the RG functions (3 g and r y m (Fig. [12]): 

l l (UR.o)=2ug R [l-^ \ . (5.21) 

9 R/ 
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9r 

lm(9R,a) = 1 -u + 2u— . (5.22) 

9*r 

Thus the theory does have a nontrivial ultraviolet fixed line = g* R , on which the 
mass operator of fermion acquires a large anomalous dimension 7 m = 1 + uj . It should 
also be noted that the anomalous dimension 7 m is continuous across the critical line 
9R = 9*R-> m contrast to the earlier phase- dependent calculation through the fermion 
propagator based on the solution of the SD gap equation. This is actually in accord with 



the suggestion p3| that the renormalization through the four-fermion Green function 
(auxiliary field propagator) in the gauged NJL model may lead to the large anomalous 
dimension in the symmetric phase as well as in the S^SB phase. 

At a = a c (uj = 0) Eqs.( |5.2iT - |5.22| ) would yield f3 g = and j m = 1. This may be 
an artifact of the symmetric renormalization based on the effective potential Eq. ( |2.67| ) 



which is no longer valid at a = a c . A possible modification at a = a c will be given in 
Section 7. 

We thus have found that the gauged NJL model is renormalized within the ladder 
approximation for non- vanishing gauge coupling a > 0. It is well known, however, that 
the pure (non-gauged) NJL model cannot be renormalized in D = 4 due to uncontrol- 
lable logarithmic divergence even in the leading approximation of 1/N expansion (or 
chain approximation). We here discuss how the existence of gauge interaction improves 
the structure of divergence. The anomalous dimension of composite operator ipif) at 
critical point of the pure NJL model is given by j m = 2. Such a large anomalous 
dimension makes the scaling dimension of ipip very small: dim ijjijj = 3 — 7 m = 1, and 
hence higher dimensional operators such as {ipip) 4 an d d^ip^d^^ip) become marginal 
operators. As a result, logarithmic divergence associated with these operators appear. 
On the other hand, the gauge interaction makes the scaling behavior of ipip softer: 
dim ipip = 3 — 7 m = 2 — uj > 1. Thus, the above higher dimensional operators becomes 
irrelevant operators and there are no uncontrollable divergence. 

In the diagrammatic picture, the softness of the scaling dimension of ipip corresponds 
to the softness of high energy behavior of Yukawa-type vertex function Eq. ( [5.14! ). 

( Fi) 

Thanks to such a soft behavior of T s , the structure of divergence is improved and the 
logarithmic divergence disappears in the presence of gauge interaction. 
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6 Operator Product Expansion 



Now that we have obtained a renormalized theory having a nontrivial UV fixed line 
with a very large anomalous dimension, 7 m = 1 + uj (gR = g R ), we can explicitly 
construct OPE both in the symmetric and the S^SB phases and see how such a large 
anomalous dimension fits in the general framework of OPE 

The OPE relevant to the fermion mass function takes the form 

- iTTT i/}(x)i)(0) = cj(p;g R ,a,m R ,ii)l 

+c^(p; 9r, a, m R , /i) [(^) R + 75(^75^] + • • • • (6.1) 

We explicitly calculate the Wilson coefficients c\ and on the nonperturbative so- 
lution which we know already. Taking the vacuum expectation value of Eq. (|6.1| ), we 
write 

-iS{p) = c t + c^{{i>il)) R ) + ---. (6.2) 
The fermion propagator in LHS takes the form 

-iS{p) = i, + ^p- + ---. (6.3) 



Comparing Eq. ( |6.3j ) with Eq.( |6.2| ), we write 

S(-p 2 ) = p 2 m R (fi)c' 1 (p; g R: a, //) + p 2 ((^) R )c^(p; g R , a, 0, //) + ••• , (6.4) 

where we have expanded the Wilson coefficients around the chiral symmetric limit 
m R = 0: 

ci(p; 9r, m Ri A*) = s + ™>rW<HW 9r^ a , y) H , (6.5a) 

pZ 

c^(P! 9r, «> m R, y) = c^(PI 9n, 0, //) H . (6.5b) 

We denote by S e xpiicit( — P 2 ) the part of the fermion mass function owing to the 
explicit chiral symmetry breaking m R ^ 0: 

S exp iicit(-P 2 ) = p 2 m ii (/x)ci(p; g R , a, fj) H , (6.6) 
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which actually takes the form 



^explicit ( ~P 



where cr exp i icit is defined by 



T s (-p 2 ; M d )(j ex pii C it, (SxSB vacuum) 

r s(~P 2 ) o"cx P iicit, (symmetric vacuum) 



(6.7) 



da. 



sol 



^explicit — C"sol ^spont 



dmr 



m H , 



(6.8) 



rrto=0 



with cr so i being the solution of the stationary condition of the effective potential Eq.( |2.62|) , 
dV(a, it = 0)/da = 0. For < a < a c it is straightforward to calculate o" ex piicit from 
Eq.( |2.73D as a Taylor series in the fermion bare mass m : 



2 — uj mo 



^"explicit 



-g (S^SB vacuum) 



2uj JL-i 
9 

m 



1 - 



(symmetric vacuum) 



9 



By comparing Eq.(|6.6|) with Eq.(|6.7|) and Eq. ( |6.9| ), we find 
12-w 1 



rf\-p 2 ;M d ) (S^SB vacuum) 



p 2uJ 9R_ { 



c'i(p;9r,v) = < 



9r 



1 1 

P 2 1 _9R 
9r 



-r?>(-P 2 ) 



(symmetric vacuum) 



(6.9) 



(6.10) 



with the renormalized vertex being given by Eq. ( |5.15[ ), where we have used 
Eq-(CT); 



m I fi\ g 



m R \ A_ 



9r 



(6.11) 
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Now, it is easy to show that the Wilson coefficients satisfy the RG equation: 



where 



= [V + 2 + 7 m (9R)]c' t (Kp;g R , a, /i), (6.12a) 
= [V + 4 - j m (g R )] c^(np; g R , a, 0, ft), (6.12b) 



V = K p g fa)—. (6.13) 
ok ag R 



Eq.( |6.12aD is readily solved to yield 

c'tiw, 9r, a; fj) - ^ (2+7 ™ ) c' 1 (p; a; y) (6.14) 
where 7^ = 7 m (5 , H) = 1 + u an d <?(/«) is the solution of 

K T^pi K ) = Pg(&)> with 9(k = 1) = flte(/i). (6.15) 



Eg .( |6.15D can be solved using the j3 function Eg. (|5.21| ): 



1 + (6.16) 



, _ ^ 1 - — ' 



9r 



9r 



Then Eg. (|6.10| ) implies a guite unusual situation, i.e., the Wilson coefficient d t does 
have a strong momentum dependence 

c'itp;^),//)-^. (6.17) 

This is combined with Eg.( |6.9| ), yielding finally the high energy behavior of £ e xpiicit: 

£ exp iicit(-«V) ^ mMp^-^c'^p; g( K ), a; //) ~ k^ 1 "^. (6.18) 

Thus the Wilson coefficient for the explicit chiral symmetry breaking term has 
a nontrivial factor Eg. fl6.10|) , yielding additional momentum dependence k 2uj , which 
actually compensates the momentum dependence arising from the anoma- 

lous dimension. As a result we obtain behavior of the fermion mass function 
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S exp iicit(— K 2 p 2 ), in accord with the solutions of the SD equation (see Appendix C): 



Vdt(-«V) - « _(1_a;) S ex piicit(-p 2 )- (6.19) 

This peculiar phenomenon with the Wilson coefficient is precisely the same as that 
in NJL£) <4 which was discovered by Kikukawa and Yamawaki [p3|l ■ Explicit calculation 
in NJL£)<4 is given in Appendix |B| where the Wilson coefficient c' t takes the form 



c't{p; git, 



1 



1 1 



D - 2 i* - 1 p 2 

9*r 

1 1 

9Y 



(S^SB vacuum) 
(symmetric vacuum) 



(6.20) 



which yields c'jQo; g(n), fi) — Solving RG equation for d t) we obtain 

c'i(«p; 9r, h) ^ c' t (p; g(n), /j,)^^^ . (6.21) 
Then in NJL£><4 we obtain 



S CX piicit(-K p) - k P m R {ii)c! t {Kp;g R ,jJL) 

~ p 2 m R (ii)c l (p] q(k), /i)K" 7 " 
= constant, 



(6.22) 



which is actually in accord with the explicit solution of the SD gap equation given in 
Appendix A. 

Let us next turn to the Wilson coefficient of (V'VOkj c^. Such a coefficient function 
can be determined from the fermion four-point function by taking x — > limit: 



i(Til>(x)$(0)1>(y)$(z) > 



I connected 



I connected 



+ 



(6.23) 

where cj does not appear in the RHS of Eq. (|6.23| ), since it contributes only to the 
disconnected diagrams. 

In the following calculation x,y and z are Fourier transformed to p,q and k, re- 
spectively. It is sufficient to evaluate the case of q = k to obtain c^. Hereafter, the 
external legs for q and k are understood to be amputated. 



42 



We first calculate the OPE coefficient function on the symmetric vacuum. The 
cx-exchange diagram (Fig. |lO|a) is given by: 



S(P^T S ( -p 2 )S(p)D aa (0)r s (-q 2 ) = -T s (-p 2 ) - 

P ~V"{a sol ) 



-T s (-q 2 ) + ■ 



(6.24) 



where we have used a relation D arT (0) = —i/V"(a so \). In addition to the above there 
exist "pure ladder" diagrams (Fig. [TTj[b) contributing to the LHS of Eq.( |6.23| ), which 
are hard to be calculated. Here we assume that such diagrams have softer high energy 
behavior than that of Eq.( |6.24 ) and ignore them in the following calculations. 

The RHS of Eq.( |6.23| ) can be calculated by (Fig. |TT|). Since a bubble diagram is 
given by the second derivative of V( qu ) = V — V( c n , we obtain (Fig. [11]) 



Z m c^(p; g R , m R = 0, n) (l - iV^ u) (a so i)D aa (0)^ T s (-q 2 
= -c^(p; 9r, m R = 0, /i)V ( / c / 1) (a so i) 



-V"(a sol ) 



Ts(-q 2 



where we have used the renormalization of the composite operator 



(6.25) 



(6.26) 



Equating Eq. fl6.24|) and Eq. ( |6.25| ), we finally obtain the OPE coefficient function 



cm(p;9r,v) 



Ts(-p 2 ) Z~ x 
P 2 Vfafa) 

r { s R \-P 2 ) G R 

p 2 N ' 



(6.27) 



where G R is defined as 



Or 

G R = A^— (6.28) 
fx 2 



We can easily see that Eq. (|6.27|) does satisfy the RG equation Eq . ( |6 . 1 2b|) . Thus 

c^(«p; 9r, p) - c$4,(p; g{t),n)K~ (4 ~ rm) - (6.29) 
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As for the S^SB vacuum, the OPE coefficient can be evaluated by replacing 
the Yukawa-type vertex in the symmetric vacuum in Eq. ( p. 27 ) with that in the broken 
vacuum: 

Tf\-p 2 -M d ) G R 

c^ip; 9r, fj) = 5 • (6-30) 

p 2 N 

It is easy to see from Eq.( |6.29| ) that the coefficient c^, having no extra momentum de- 
pendence, yields a correct high energy behavior of S dyn (— p 2 ) = p 2 {(4>ip)R)c 1 p 1 p(p; gn, a, //): 

E dyn (-K 2 p 2 ) = KY((^) R )c^(p;g( K ),fi) ~ k~^) „ K -d-») ) (6.31) 

which indeed agrees with the solution of the SD equation given in Appendix C. 



This result is also similar to that of |[23|| . See Appendix B for details 



7 M-Dependent Renormalization 

In the previous section, we have carried out the symmetric renormalization of the 
gauged NJL model with a < a c (u 7^ 0), based on the effective potential around the 
symmetric vacuum Eq. (|2.67D . Such a renormalization cannot directly apply to a = a c 
(lo = 0), since the expansion Eq. (|2.67|) loses its validity at this point. However, break- 
down of the expansion Eq. fl2.67p does not necessarily imply the non-renormalizability 
of the gauged NJL model at a = a c . Actually, as we have studied in section 3, the 
fermion scattering amplitudes remain finite even at a = a c as well as at a < a c , once 
the bare four-fermion coupling is fine-tuned so as to make the fermion mass finite. 
This fact suggests that the renormalization may be possible by a suitable definition of 
renormalized parameters even at a = a c . Actually, the general form of the effective 
potential Eq.fl2.62j) remains valid even at a = a c . 



In this section we present yet another renormalization scheme which renormalizes 
the effective potential Eq.( p.62[) instead of Eq.( 2.67 ). The renormalization described in 



this section introduces a redundant mass parameter M (M-dependent renormalization) 
and M = corresponds to the symmetric renormalization. The M-dependent renor- 
malization is related to the symmetric renormalization via a finite renormalization at 
< a < a c and remains valid at a = a c . 

The Yukawa-type vertex Eq. ( p. 85 ) at the non-trivial M 7^ is renormalized by the 
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wave function renormalization: 



1-w 



V » J \ v , 



l+UJ ' 



(7.1) 



which leads to the renormalized Yukawa-type vertex in the S^SB vacuum M = M ( 

dX{-p 2 )/dM\ 



Tf\-p 2 ,M d 



I M=M d 



(2-w) 



Ci 



M 
. . 



1-w 



£>1 



M 



M 



l-w 



(7.2) 



Note that the same renormalization condition simultaneously makes finite the Yukawa- 
type vertex in the symmetric vacuum: 



rf {-p 2 )=Tf (-p 2 ;M = 0) 



-P 



2N -(l-w)/2 



1 + UJ V ^ 2 



Ci UJ 



(7.3) 



We next define the renormalized four-fermion coupling and renormalized fermion 
mass rriji by: 



and 

mo 




«(?) 

+Dl A. 



fDl (i, 
(t) 



// 



V / 



// 



\9 9 J 



1_1 



M 9 R , 



-(!-«)■ 



(7.4) 

Here g*^ and c/j=, are arbitrary parameters corresponding to the choice of M in this 
renormalization scheme. It is easy to see that the Eq. ([7.4Q reduces to the symmetric 
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renormalization Eg. ( |5.17| ) in the M — > limit. Also note that unlike the symmetric 



renormalization, in Eg .(7.4) does not corresponds to a "critical" coupling, unless 
M/n = 0. 

Starting from the bare effective potential Eq. (|2.62j ), we obtain a renormalized ef- 
fective potential in the A — > oo limit: 

^ 2 V R (a R ,7r R = 0) 1 m R a R 1 / 1 1 \ f aR \ 2 



N 



9r V 



1 I 1 

+ 2 r 



2 V^s 05/ V 



R xr, 



'M\ 2w 2-u 



Ci 



Ci 



. , 



G* 













2 - LU 






2 C ' D ' 











Ci 



Eq.( |2.40| ) and Eq.( [7.1| ) lead to the relation between M and the field a R : 

C 1 + D 1 



M 



2-oj 



(7.5) 



Note that M is a certain constant but not a field like M. 
We choose the parameters g*^ and g^; 

1 1 

= h regular function of w, 

9*R -9R u 



(7.6) 



so as to take a sensible uo — > (a — > a c ) limit. Hereafter we take the simplest choice: 



9r = 9 



:i-u) 



(7.7) 
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Such a choice of g* R and g* R leads to the following effective potential at uj = (a = a c ): 
4n 2 V R (a R ,n R = 0) 



N 



2 3 uM 



1 + S - In 



M 



2 ' 



(7.8) 



The stationary condition of the effective potential Eq. ( |7. 5| ) yields the (renormalized) 
gap equation. Solving the gap equation in the chiral symmetric limit m R = 0, we obtain 
the scaling relation of the dynamical mass of the fermion = (M): 



M£ - 
M 



D t ( M' 



2^ 



1 1 








2 


(?) 






9s 9R 









1 1 

~9*R 9* R 



The critical NJL coupling g c £ l is determined from Eq. ( [7.9| ) at = 0: 

2w 



1 



1 



o- rit a* 

9r iJR 



+ 



M 9r, 



2 ■ 



Note here that the critical coupling can be calculated as 



9t 



g% forM/// = 
g* R , for M/jj, — > oo 



(7.9) 



(7.10) 



(7.11) 



If we take M = (a natural choice in the S^SB phase) in Eq. (|7.9|), we find: 



Md 



2^' 



1 1 

'A J__ J_' 



(7.12) 



Now, we are ready to discuss the RG properties. It is straightforward to derive the 
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(3 function from the definition of renormalized four-fermion coupling Eq.( |7. 1| ): 



Pg(j9R,a;M/fi) = -2ug R 



i— L )c! (H 

2 \9j SrJ V A , 



-2ui 



1 M d? /m 



2^ 



^ 9R, 



M 



and the anomalous dimension 7 m from Eq.( [7.4|) : 



Hm(9R,oi;M/fi) = 1 + w- 



-2a; 



-I -Wff 



(7.13) 



(7.14) 



Putting M = in Eq.(fTT3D and Eq.(|7TT4l), we find (Fig. 0) 



&(^a;M = 0) = 2^ ^1 - p 

7m (^,a;M = 0) = 1-^ + 2^, 

2j? 



(7.15) 
(7.16) 



which coincide with Eqs.( [5.2l| - |5.22| ) as they should. 

Plugging Eq. (|7.12|) into Eq.( [7.13j ) and Eq. (|7.14|) , on the other hand, we obtain 
(Fig. 



P g (g R ,a;M = M d ) 
lm {g k , a;M = M d ) 



- 2 (9r~ 9*r){9r- 9*r), 



a 

2g R + 77-, 

2a r 



(7.17) 
(7.18) 



which remain valid even at a = a c and coincide with the form of the RG functions of 
bare parameters, Eq.( |U.18| ) and Eq.( |C.23| ), given in Appendix C. 

Finally, we make a brief comment on the relation of the symmetric renormalization 
and the M-dependent renormalization. From Eq. (|5.13|) and Eq. (|7.1|) , we obtain 



VR = Z a a R ; Z a l = l 



D l ( M 



2uj 



(7.19) 



in the A — > oo limit. Using Eq. fl5.17| ) and Eq.(T^), we obtain a relation between the 
renormalized four-fermion coupling of the symmetric renormalization and that of the 
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M-dependent one: 



9* 9R ° \9* 9r) Ct\^J \g* g* 



where we have used g* = g* R = g* R and g* = g%. It is easy to show that the current 
mass in the M-dependent renormalization scheme is expressed by 

-_,gn 

m R = Z m m R , Z m = Z a — (7.21) 

9r 

in the A — > oo limit. Thus, the symmetric renormalization and the M-dependent 
renormalization are connected by a finite renormalization unless uo = 0. 

The OPE coefficient functions in M-dependent renormalization c' t and are cal- 
culated from 

c't(p; 9r, a,M/fi;n) = Z m c' t (p] g R , a; /i), (7.22a) 

c^(p; gR, a, M/n; //) = Z^c^p; g R , a; //). (7.22b) 



By using Eq. (glQp , Eq. (lQ7|) and Eq. flOOl) , we find 



12 -"rf(V;M d ) 



c'i(g R , a, M/n\ n) = — P ^ D i / / x Y\ 7 ^' 23 ^ 

c^(g R ,a,M/^^ = -l^Mrf\-p 2 ;M d ) (7.23b) 



for the S^SB vacuum, and 



irf(V;M = o) 

c'i(gR, a, M/fi; (j) = £ ^ (7.24a) 

l_^_Z 2 ^i| — ] fJ__J_] 

9* a Cl\ii) \g* g*J 

c^(g R ,a,M/fi ]f i) = -^^^Tf\-p 2 -M = 0) (7.24b) 
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for the symmetric vacuum, where we have used 



rf \-p 2 ; M) = Z-'Tfi-p*- M). (7.25) 

We notice that 

Z (j -> i (7.26) 

-2cu\nf + 0(lu 2 ) y ! 

as u — > 0. Therefore on S^SB vacuum and c'j remain finite in u — > limit, while 
in the symmetric vacuum diverges. 

Special attention should be paid to the uj — > limit where two renormalization 
schemes are not connected by a finite renormalization. The above calculation implies 
that the OPE coefficient function cannot be made finite in a vacuum-independent 
manner. Such a peculiar phenomenon may be an artifact coming from the infrared 
singularity in the symmetric vacuum at a = a c . 



8 Conclusion and Discussion 

We have presented renormalization of the simplest gauged NJL model, gauge theories 
with standing gauge coupling plus four-fermion interaction, in the ladder approxima- 
tion. Through the CJT effective potential written in terms of the auxiliary fields, we 
have established, in an analogous manner to NJLd<4, a phase-independent renormal- 
ization (symmetric renormalization) valid both for the symmetric and the S^SB phases 
as far as a ^ a c . Accordingly, the (3 function and the anomalous dimension were ob- 
tained phase-independently in both phases: The theory has a nontrivial UV fixed line 
and a large anomalous dimension there. The OPE was explicitly constructed, which is 
consistent with the large anomalous dimension in both phases. The Wilson coefficient 
for the unit operator has an extra power behavior other than the anomalous dimen- 
sion. The symmetric renormalization done on the symmetric vacuum breaks down at 
the end point a = a c of the critical line, while the M-dependent renormalization still 
remains valid. 

The reason why the renormalization is possible in NJL£) <4 and gauged NJL model 
is very simple. In NJL£> <4 the four-fermion operators become relevant /marginal: 
dim ((?/5-?/>) 2 ) = 2dim('0'0) = 2(D — 1 — 7 m ) = 2 < D, while other fermion oper- 
ators become irrelevant: dim ((ijjij)) 4 ^ = dim ^9 /1 ('0'0)9 M ('0'0)^ = 4 > D, etc. for 
2 < D < 4. In the same way, the four-fermion operators in the gauged NJL model 
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become relevant /marginal: dim {[^)il)) 2 ^j = 2(3 — j m ) = 2(2 — yjl — a/a c ) < 4, 



while other operators become irrelevant: dim ( (ipip) 4 ] = 4(2 — Jl — a/a c ) > 4, 



dim (d^VO^W)) = 2 ( 3 - V 1 ~ «/«c) > 4, etc. for a c > a > 0.0 

At this point it should be emphasized that our renormalization breaks down at the 
pure NJL limit a — > 0. In that limit our renormalized effective potential Eq. ( |5.19| ) has 
a singularity of l/a,0 a signal of the appearance of the logarithmic divergence for the 
four-point vertex of the local auxiliary fields. Such a divergence can only be removed 
by introduction of eight-fermion operator ( "A0 4 counter term" ) which now becomes a 
marginal operator; dim ((?/>?/>) 4 ) = 4(2 — y 1 — ot/ct c ) — > 4 at a — > 0. The 1/a sin- 
gularity also appears in the auxiliary field propagator Eq.( p.2(J| ), another signal of the 
logarithmic divergence in the induced kinetic term of the auxiliary field. This again 
can only be removed by introduction of derivative-type four-fermion operator ( "counter 
term for kinetic term d^(f)d IJ, (f) n ) which also becomes marginal; dim (^(^/^c^ (?/>?/>)) = 
2(3 — Jl — a/ot.c) —* 4 at a — > 0. Without such extra "higher dimensional" operators, 
the pure NJL model in four dimensions cannot be renormalized even in the nonpertur- 
bative sense of 1/N expansion. The presence of gauge interaction turns these "higher 
dimensional" operators into irrelevant ones and hence make the renormalization pos- 
sible without such additional "counter terms". [^2], |15|, [16| Similar comments also apply 
to the NJL/X4: The above 1/a singularity corresponds to the 1/{D — 4) singularity in 
Eq.(p|) and Eq. fl5T2| of NJL D<4 . 



Such a formal resemblance between the renormalization of NJL£> <4 and that of the 
gauged NJL model may not be a mere accident. As is well known, the renormalizable 
NJL£><4 model is equivalent to the Yukawa model with the Yukawa coupling lying on 
the nontrivial IR fixed point in the 1/N leading order. Actually, the j3 function of 
dimensionless Yukawa coupling y is given by 



,D \ r(3 — D/2) (T(D/2)) 

My )= (il_2)y + 4JV— ' y 3 , (8.1) 

1 2 J (4n) D / 2 T(D - 1) 



8 These higher dimensional operators can be explicitly shown to be suppressed in part by the power 
of the cutoff in the low energy physics through renormalization. 

9 The 1 1 a singularity instead of the logarithmic divergence is an artifact of the inadequate approx- 
imation taking only the first two dominant terms in the solution of SD equation for a > 0. Taking 
account of the third dominant term, we can recover the correct NJL limit |14|, [2q|. 
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which indeed has a nontrivial infrared (IR) fixed point 



2 w /2 np - 1) 

V = tt^t, , ; ^r, ^ ( 8 - 2 ) 



4NT(2-D/2) (r(D/2)) 

corresponding to the NJL£> <4 model. Although two parameters (Yukawa coupling and 
scalar boson mass) are needed in order for the Yukawa model to be renormalizable, the 
number of parameters is now reduced by this constraint to that of NJL£> <4 . This is the 
essence of the renormalizability of NJLd<4- This argument obviously breaks down at 
D = 4 where the nontrivial IR fixed point disappears. This situation reflects the fact 
that the NJL model is not renormalizable in four dimensions as an interacting field 
theory even in 1/N expansion. 

However, in the presence of gauge interaction (a ^ 0) the Yukawa model shows a 
similar structure to Eq. flS.lfl even in four dimensions: 



6C F a N 

— y + — 2 



p y (a, y) = —y + —y\ (8.3) 



where Cp is the quadratic Casimir of the fermion representation and we have ignored 
the graphs containing scalar particle loop. Now, the (3 function of gauge coupling may 
be parameterized as 

"iCpjn 

(3 a (a } y) = a% (8.4) 

A 

where A(= 18Cp/(llN — 2Nf) for iVy-flavored SU(N) gauge theory) is the measure of 
the running speed of the gauge coupling, i.e., large A(^> 1) means "walking" (slowly 
running) gauge coupling and A — > oo corresponds to the non-running ("standing") 
case studied in this paper. For A > 1 the above (3 functions lead to the RG invariant 
IR stable subspace of the couplings: 

y = (8.5) 



in analogue to Eq . fl8.2p . This is another expression of a suggestion fig, |4| that the 
presence of gauge interactions with A > 1 may make the theory renormalizable. The 
explicit renormalization procedure in the present paper is actually a concrete example 
of this suggestion in the special case A — > oo. Such a possibility in the case with 
running/walking gauge coupling A < oo will be further studied in the subsequent 
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paper 
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A NJL^<4 model 

In this appendix, we derive the effective potential and the auxiliary field propagators 
in NJLd<4 model. 

The auxiliary field technique enables us to rewrite the original lagrangian Eq. fl5.ip 
into an equivalent one: 

C = - ip(a + ni^ip - V^a, vr), (A.l) 

where classical potential of auxiliary fields a, ir is given by 

N 



V (d) = q 



7T 2 l — mod 



(A.2) 



The effective action is evaluated solely from the fermion integral within the 1/N leading 
approximation: 

r[er, 7r] = — zlndet (0 — a — 27571- ) — J d D xV( c i) (c, tt) . (A. 3) 

The effective potential V is calculated from the effective action 

rjer = const, 7r = const 1 
V(a, vr) = - J - 1, (A.4) 

where Q is the space-time volume. Thus, we obtain 

V(a, tt) = V {cl) (a, vr) + V (qu) (a, tt), (A.5) 

where 

rA d D p ( a 2 + vr 2 ' 



/•a a p / a -+- 7r \ 

jfdpl(pir*-Hn(l+ a -^Y (A.7) 



(47r) D/2 r(D/2) 

with the momentum integral being regularized by the ultraviolet cutoff A. Using the 
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following formula 



dtt 



z-l 



In 1 + 



— B(z -1,2- z) 



(A.8) 



it is rather straightforward to integrate Eq.( [A.7| ): 

(47r) D/2 rp/2), 



4N 



-V(qu)(^7r) 



A D—2 2,2 /- ( '2 , 



7T 



D/2 



9 



2 - D/2 



D 



E 



A 



D i 



<7 + 7T 



n+2 



o(D-4-2n)(n + 2) V A 



(A.9) 



with = D/2 - 1 and ( D = B(D/2 - 1,3 - D/2). Finally, we find the effective 
potential 



(Au) Dl2 Y(D/2) 



4NA D 
1 m a 



V(a,7r) 



1 1 \ a 2 + n 2 



, 9 0\ D / 2 

1 CzW ff +^ 



5 A 2 



JJ 



9 2A< 



D 

2 



D 



iV 



+ o 



' 2 i 2 \ 2^ 

(7 + 7T ) 



A 2 



,(A.10) 



with the dimensionless four-fermion coupling g being defined by 

4A D - 2 



9 = 



(4n) D/2 T(D/2) 



G. 



(A.11) 



It is easy to see that the last term ~ 0{A D ~ i ) in Eq. flA.lO ) becomes negligible for 
sufficiently large A and will be neglected in the following calculation of the propagators 
of auxiliary fields. Note that g* — > and (d —> oo as D — > 2, so that the divergences in 
the second and third term cancel each other to give a well-known logarithmic potential 



in the Gross- Neveu model p7f: 



7T 



1 



-V(a,n) = -m a-- ^ 



la 2 + ir 2 cr 2 + 7r 2 



+ 



1 - In 



V + TT 2 ^ 

A 2 



+ A 2 0{ 



' 2 i 2 \ 2 



A 2 



(A.12) 
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Let us now consider solution of the gap equation dV/da = for Eq.( [A.10l) : 



lm / 1 1\ a 
°"^T + W~ ~g) A 



C 



D / (7 



2 v 

2 



D-l 



(A.13) 



It is convenient to parameterize the solution of the gap equation a so i as: 

"sol ""spoilt "explicit; 

with cr spont being the solution in the chiral limit (m = 0): 

1 1 



(A.14) 



Cd 



2 - D/2 



-a. 



D-2 
spont 



A 



D-2 



.9 9. 




for S^SB vacuum (g > g*) 
for symmetric vacuum 



(A.15) 



Actually it is easy to see from Eq. ( |A.10|) that the symmetric vacuum (a) = (ir) = 
becomes unstable in the strong coupling region g > g*. In this vacuum the fermion 
acquires the dynamical mass 



S dyn(-P 2 ) = r 5 (>, q = 0)cT spont = CT spont , 



(A.16) 



where T s (p, q) is the Yukawa- type vertex defined in an analogous way to Eq.( [2.10| ) and 
calculated as 

r s (p,q) = l. (A.17) 



The nonvanishing bare fermion mass (m ^ 0) causes non-zero cT cxp ii C i t : 



da, 



sol 



a, 



explicit "sol "spont 



dmn 



m + 



(A.18) 



mo=0 



56 



We can evaluate explicit as a Taylor series from Eg. ( |A. 13| ) in m : 

1 m 



^"explicit 



D-2JL-1 

9* 



m 

9* 



+ ■ 



for S^SB vacuum 
for symmetric vacuum 



(A.19) 



Hence the current mass of the fermion may be written as 



^explicit( — P 2 ) — ^s(p, q — 0)cr cxplicit — cr cxp i icit . 



(A.20) 



Thus, the current mass S exp i icit has the same (constant) high energy behavior as that 
of the dynamical mass in the NJLb <4 model. 

The auxiliary field propagators D~l and are given by the second derivative of 
the effective action: 



S 2 T[a, tt] 



5a(x)5cr(Q) 
5 2 T\a, tt] 



5ir(x)5ir(0) 



0"=0" S0 1 ,7T=0 



<T=(T ao l,7r=0 



(A.21a) 
(A.21b) 



They are evaluated as 



d D k 
(2ti) d i 

d D k 
(2n) D i 



-tr 



-tr 



1 % 
-i- 



ft - er so i ft~i~ cr sol 

i i 
75- 75" 



a 2 V (c i)(<j,7r) 



ft - O-sol ft - i - CTsol 



da 2 

'(d) v 



(A.22a) 



f=<T st ,i,7r=0 



d 2 V( r ])((T, TT) 



dn 2 



(A.22b) 



C=<T sol ,7T=0 



The subtraction at zero momentum makes these expressions finite in the A — > 00 
limit: 10 

iD£(-?)-iD£{0) 



10 Thus, we do not need to worry about keeping the regularization chiral invariant for the auxiliary 
field propagator with non-vanishing momentum q 2 7^ 0. 
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r(2-£>/2) f /Y 



8A— ^— < | — - ^soi I / dx 



(4tt) d / 2 I \ 4 
T(2-D/2) (q 



^soi - x(l - x)q" 



D/2-2 



af fM (A.23a) 
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(4tt) d/2 I 4 j Jo 



— / dx 



D/2-2 



(A.23b) 



The second derivative of the effective potential gives the auxiliary field propagator at 
zero momentum: 



iD-^-q l = 0) 



d 2 
da 

d 2 



it 



tD-^-q 2 = 0) = -—V(a,n) 

07T 



°"=' T soli 7I "=0 



<7=CT so l,Tr=0 



(A. 24a) 
(A.24b) 



For the symmetric vacuum <7 so i = we find 
iD-l{-q 2 = 0) = iD-l{-q 2 = 0) = - 



4A A D - 2 1 1 1 
(47r) D / 2 rp/2) \g~~g~* 



(A.25) 



Combining Eq. (|A~23D with Eq.( [A25l) , we find 

(47r) D / 2 



4A 
(47t) d / 2 



4A 

with £d defined by 



V{D/2)iD^{-q 
T{D/2)iD-l{-q 



— I ; | _j_ | \/> 



D 



,9 9 
1 1 



A^" 2 + 



2 - D/2 

_Jd_ 

2 - D/2 



(-gT /2 -\(A.26a) 
(-g 2 ) I?/2 - 1 ,(A.26b) 



5(3 -D/2, D/2 - 1) 



D 



(A.27) 



T(D-l) 

For the S^SB vacuum a so \ ^ in the strong coupling region g > g*, the auxiliary field 
propagators at zero-momentum are evaluated as 



iD- A (V = 0) 



8A 



(4tt) d / 2 



r(2 - D/2)^- 2 , ^(-g 2 = 0) = 0, (A.28) 
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which lead to 



,DJ2 



^^r( J D/2)^ CT 1 (-g 2 ) = 2(^-^1 fdx 
AN \ 4 / Jo 



(47T 



.D/2 



AN 



-T(D/2)iD^(-q*) = 2 



Q 



4 ) Jo 



dx 



^sol - x(l - x)q 2 
D/2-2 



D/2-2 



(A.29a) 
(A.29b) 



Thus, the auxiliary field a acquires a pole at q 2 = Aa 2 ol = 4m 2 , while n becomes 
massless NG boson. Note that m a = 2m is independent of D. 



B OPE in NJL D<4 model 

Let us consider the time ordered bilocal operator T ip^x^ifi) 



. The OPE for it reads 



i)(x)ip(0) =c 1 (p\g R ,m R ,fi)l+c^i ) (p]g R ,m R ,fi) (ipip) R + TsC^TsVO 



)R 



+ ■ 



It is useful to expand c% by m R : 



ci{p; 9R: m R , n) = — + m R d t {p\ g R ; fi) + 



Corresponding to Eq.( |6.6|) , we have 



p- 



S cxp iicit(-P 2 ) = p 2 m R (fi)c' 1 (p; g R ; fi) + ■ 



From Eq. flOij ) and Eq. (pCT9| ) we obtain 



1 Z m m R 

-g 1 + • ■ ■ for S^SB vacuum 



S eX plicit( — P 2 ) — i 



D-22--1 

g 

ZmtTlR 



g* 



+ ■ 



for symmetric vacuum 



where Z m = m /m R is given by Eq.(pT 



D-2 



m _ g I fi 
m R g R I A 



(B.1) 



(B-2) 



(B-3) 



(B-4) 



(B.5) 
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Eq.(|EjjP and Eq.QET§) lead to 



1 



9 



1 



9r 



(B.6) 



F 9r 

Comparing Eq. (|B.3|) with Eq.( B.4j ) and Eq. (|B.6|) , we obtain 11 

1 1 1 



9*R 

1 1 

9*R 



for S^SB vacuum 



for symmetric vacuum 



(B.7) 



Note that the Wilson coefficient d t in the S^SB vacuum remains finite in two di- 
mensions, while c'j calculated in the symmetric vacuum vanishes in D — > 2. This is 
consistent with the one phase structure in two dimensions. 

Let us next consider the coefficient function of i^if), c^. This can be calculated 
from the fermion four-point function by taking x — > limit: 



-i(Tty(x)$(0)i/>(y)${z))) 



connected 



connected 



The Wilson coefficients of the unit operator do not appear in the RHS of Eq. ([B.8|) , 
since they only contribute to disconnected diagrams. 

In the following calculation x,y,z are Fourier transformed to p,q,k, respectively. 
The calculation at q = k is sufficient to determine c^. The LHS of Eq.( B.8|) is 
evaluated: 



S{p)iT s (- P 2 )S(p) 



1 



-V"(a sol ) 



*r s (-q 2 ) 



1 



1 



P 2 -V"(a sol ) 



r s (-q 2 ) + 



(B.9) 



where the legs for q, k are amputated and the scalar vertex is defined by T s 
The RHS of Eq. ( B.8|) can be evaluated by Fig.[y] which reads 

Z m cj,Jp;g R ,m Rj fi) fl -iV/'Ja^D^O)] T s (-q 2 ) 



1. 



ii 



This agrees with Ref. [|23| except for the factor !/(£) — 2) in the SxSB vacuum. 
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= -Cfa (p; gfl, m fl , /i) V ( ci) (o-soi) ~ m — -r 5 (-g 2 ), (B.10) 

-»/ (O"sol) 

where in the last line we have used Eq.( |A.241 ). 

By comparing Eq.( |B.9| ) and Eq.( [B.10|) , we obtain the Qw, |f23 



c^a(p; ^Wij = 0,/i) 



7-1 



(B.11) 



P 2 ^(^sol) 

1 (47r) g / 2 r(D/2) ^ 
p 2 4iV /i D ~ 2 

p 2 N ' 

where G# is defined by Gr = (4:Ti) D ^ 2 T(D /2)gnfx 2 ~ D /4 and we have used 

J_ = G = (47r)^ 2 r(ZV2) g 

C Renormalization Group Functions from the So- 
lution of Ladder SD Gap Equation 

Let us consider the RG of bare parameters which was first discussed by Miransky || in 
the analysis of the ladder SD equation of dynamical mass of fermion in the quenched 
QED. In this argument the bare parameters of the theory are required to depend on 
the cutoff A so as to fix the fermion dynamical mass M<j, with such a flow of bare 
parameters being identified as the RG evolution. 

Solving the ladder SD gap equation of the gauged NJL model, Kondo, Mino and 
Yamawaki [|BJ and independently Appelquist, Soldate, Takeuchi and Wijewardhana 
13] obtained the critical line 



9 = 9 



\(l + cu) 2 , u = y/l- a/a c , (C.l) 
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and the scaling relation near the critical line: 
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/ 1 1 \ 1/2uJ 
Cl £ g 

V 9* 9 J 
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cxp 



i r 

9* 9 
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(J/2 



5- 



g-{l-^)/A< 



U! 



(0 < a < a c 



(a = a c ) 



{a > a r 



(C.2) 



with n = 1 being the ground state solution, where g* is defined by g* = (1 — co>) 2 /4. 



The usual SD equation |]TT| is obtained from the stationary condition of V[E, E 5 ] 



5 



5E(p : 
5 



*£5(l£) 



f[£,£ 5 ], 



(C.3) 
(C.4) 



which read 



E(p 2 ) = m + — J dp 



A 2 



_ _L_ / ^fc _ 

B p 2 +E 2 + E 2 E fc 2 + E 2 + E?~ 



tf(z4A£)(C.5a) 



p 2 £ 5 Cp 



p 2 + E + Ej 



+ 



dfc 2 



^ B E5(/c f 



fci + S 



K(p 2 ,A; 2 ). (C.5b) 



Since E 5 can be rotated away by the chiral symmetry, it is sufficient to consider the 
SD equation 



E(p|) = m + J dp 



A2 , ,vlWi 



h 



a2 „ a gg(g; 

fc* + E 2 



dki 



K(p 2 E ,k 



2 'i), 



(C.6) 
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which is equivalent to the differential equation 



d 

dpi 



+ 2 



d 3C F ol 
dp 2 E 4vr pi + E 2 (p 



E(p : 



(C.7) 



with IR BC: 



and UV BC: 



d 

i im p e —t e (Pe) = 0, 

pI^o dpi 



1+1 + 



3C F a/AnJ^ E dp 



m . 



(C. 



(C.9) 



Irrespective of presence or absence of the explicit fermion mass mo, the solution of this 
differential equation exhibits the same asymptotic behavior in the high energy region: 



M 



c-i 



M 2 



-(l-w)/2 



+ d 1 



P^ 
M 2 



-(l+w)/2 / o \ -3(l-w)/2-l 



(CIO) 

where this time M is not a field but a constant having the dimension of mass and 
gives the mass scale of the solution. For details on the behavior of this solution, see 
section 2. 

First of all, we calculate the mass renormalization constant Z m d la Miransky ||. 
From the UV BC, the equation of state is obtained: 



_ /My 

a va; 



9-9* n (M 



g* V A J 



(en; 



Then we obtain 
~dM ~ 



' g* V A 



(C.12) 



In the chiral limit m = 0, the dynamical fermion mass obeys the scaling law 
Eq.(CJ) 



Ma 
A 



~g*c l9 -g* ~ 

g*D 1 g-g* 



i 

2jJ 



(0 < a < a c ). 



(C.13) 
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Then the renormalization constant in the region < a < a c is obtained as 



7 = dm ° \ 

~~ <9M ' M=M <i 



Thus we obtain the mass renormalization constant 



2— k\ (g - g*){g - g*), 



\ A \l{9-~9*){9-9*) M d 



Z m — < 



A 



(0 < a < a c ) 
(a = a c ) 



2A 



,\J( 9-9 c )(9-9 c ) M d 
Vl + u' 2 A 



[a > a r 



where g c = (1 + icu') 2 /4 and g c = (1 — icu') 2 /4. 

Since the bare chiral condensation (?/>■?/>) is given by 
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4vr 2 7o ™ p l + £ Vb ) 

~47rVl-^[G?-r)(£-S*)] 
iVA AM j 



1/2 



(0 < a < a c ) 



4tt 2 



r 

NA' 



AM 



[a = a r 
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{ 47T 2 VTT^[(g-g c )(g-g c )] 1/2 
the renormalized condensate {(^iP)r) near the critical line is calculated as 

NA 2 



{(W) R ) = Z m {W) = 



NAl , 



2tt 
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(a = a c ) 



M| (a > a c ) 
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This shows that ((■0 , 0)_r) is g-independent and the RG flow can be identified with 
the fixed-a line (upward direction) |TJ, |TJ, [OJ . 

Once the RG flow is so identified, the j3 function of bare four-fermion coupling was 



explicitly calculated from Eq. ( |C2|) in the gauged NJL model pO, |21 



do 

P g (g, a) = A^-L, Md = -% " 9*)(9 ~ 9*), (9 > 9*)- 



(C.18) 



The anomalous dimension is obtained [[KJ 



lm{g,oc) = -A 



d In Z m 

dA a ' Md 



-A 



d In m , 
<9A 



1 + u, (g = g*). 



(C.19) 



Now we calculate the anomalous dimension above the critical line g > g* from the 
solution of the SD equation. The scaling law leads to 



9-9 = 
with F being defined by 



1 F 

Y^p^*-g), 9-9* = y^f^ 9 ~~ 9 >> 



g* Ci V A 



(C.20) 



Then we get 
d\n(g-g*) 



d In A - M d d{ A 



A dln(l-F) F 



1-F 



-^ l r J ^ = -^9-g*)- (C.21) 
g -g 



Similarly, we get 
d\n(g-g*) A d[\nF — ln(l — F)] 



din A 



M, 



A - 



-2a; - 2a;- 



1 — F 



-2w-2(p-/). (C.22) 



Accordingly, we obtain the anomalous dimension: 

dlnZ m 



din A 



i ~±J^Mg-g*)+Hg-g*) 



2<91nA' 



a 



l+u + 2(g-g*) = 2g + . 

2a r 



(C.23) 



This is shown to be valid also in the region a > a c and coincides with the anomalous 
dimension numerically obtained in the gauged NJL-model with running gauge coupling 
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m. 



D Renormalization Group of Bare Parameters 

In the previous appendix, we have defined and calculated the RG function of bare 
parameters. It is evident, however, that such a definition of RG function is applicable 
only in the SxSB vacuum of strong coupling phase g > g* where fermion acquires 
non- vanishing dynamical mass. 

Thus, we define here the RG of the bare parameters by using the effective potential 
expressed in terms of mass scale of the fermion M. This enables us to obtain the RG 
functions also in the symmetric phase. 

We first consider the case of mo = 0. The effect of explicit chiral symmetry breaking 
term will be discussed later. In section 2, we have obtained the effective potential: 



•>7T 



V(M) 



N A 4 



1 



,9 9, 



Ci 



Ci 



M 
A 



4-2c; 





\ 4+2w / \ 4' 

M\ 2 + uj ( M 

D 1 | — + G x — 

A 2 I A 



Actually, the stationary condition of Eq.( p.l| ) dV/dM = has a nontrivial solution 
M = Md 7^ for strong coupling region g > g* and leads to the correct scaling relation: 




Cl g* g 



Di 



1 ' 



(D.2) 



9* 9 



Let us define the RG flow of the bare four-fermion coupling g so as to make Eq.( p.l 
independent of A: 12 



8vr 2 dV(M) 

= — = 2u 

NA 3 dA 





4+2o;~ 



12 It is evident, however, that the effective potential Eq.( D.l ) cannot be made A-independent for all 
the region of M. Thus, we need to specify the value of M to define the RG flow of bare parameters. 
The plausible choice of M is the value of the stationary condition. In this sense, the RG flow of the 
bare parameters depends on the choice of vacuum. 
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9 2 U 



M 



M 



where (3 g is defined by f3 g (g, a) = Aj^. Namely, the j3 function reads 

1 1 



Pg(9,<*) = ~ 2uj 9 



2 \9 9, 
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,9* 9, 
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(D.3) 



(D.4) 



We first consider the symmetric vacuum M = 0. In this case, Eg. ( p.4p becomes 



Pg{g,a) = 2ug 1 



(D.5) 



which corresponds to the phase-independent RG evolution of the renormalized four- 
fermion coupling discussed in section 5. 

Another choice of M is the solution of the gap equation Eq.( |D.2| ) in the strong 
coupling region g > g*, which gives 



P g (g,a) = -2(ci-g*){g-g*). 



(D.6) 



Eq. (p.6|) agrees with the original definition of the RG flow of the bare parameters 
Eq.CT). 

These two results Eq. ( p.5|) and Eq. ( p.6|) give the same result for the region g — g* 
where M is sufficiently smaller than cutoff A. On the other hand, the deviation becomes 
significant when the four-fermion coupling g is far beyond the critical line g* and M/A 
in Eq.( p.2|) is not negligible. 

For the anomalous dimension 7 m , we need to evaluate the explicit chiral symmetry 
breaking term proportional to the bare mass of fermion mo in the effective potential: 



explicit 



(M) 



A^A 2 m a 
Att 2 g 



(D.7) 
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The result is: 



8vr 2 Explicit (M) 2m 
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N A 4 A# 

The anomalous dimension 7 m m = — A<9m /<9A can be obtained in the same way as 
the (3 function; aVe ^ licit = 0. We find 



cl\ 















-2u 
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(D.9) 



Again, the definition of the anomalous dimension 7 m depends on the choice of M. The 
symmetric vacuum M = gives the anomalous dimension 



g 

lm(g,oc) = + 2a;—, 

9* 



(D.10) 



which takes the same form as Eq.( |5.22| ), the anomalous dimension for the renormalized 
coupling in the symmetric renormalization. It should be noted that the anomalous 
dimension 7 m is continuous across the critical line g = g*. On the other hand, using 
the value of M in S^SB vacuum Eq.( D.2 ), we find 



a 



lm(g,oi) = 2g H . 

2a c 



(D.ll) 



Eq.( p.ll]) for g > g* can also be derived by use of the SD gap equation, as shown in 
Appendix 0. Eq.( p.lO|) and Eq.( p.li|) give the same anomalous dimension Eq. (|C.19 ) 
at g = g*. 
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Figure Captions 



1. Critical line in (a,g) plane. It separates the spontaneously broken phase (S^SB) 
and the unbroken phase (xSym) of the chiral symmetry. 

2. The lowest order diagram in the two-particle irreducible part k 2131 ^] of the CJT 
potential of the gauged NJL model written in terms of auxiliary field Eq.(2.2). 
The wavy line and the solid line with shaded blob represent the bare gauge boson 
propagator D^ v and the full fermion propagator S, respectively. 

3. The lowest order diagram in the two-particle irreducible part of the CJT potential 
of the gauged NJL model (without auxiliary fields), Eq.(2.1). 

4. Momentum assignment of Yukawa type vertex T s (p,q). Dashed line represents 
auxiliary field propagator. 

5. Amputated multi-point Green function at zero momentum of a. External fermion 
lines are amputated. The shaded blob stands for the induced a vertex (V^ n \a)). 

6. Auxiliary field propagator. The solid line with shaded blob represents the full 
fermion propagator. 

7. The second derivative of the auxiliary field propagator. A slash represents once 
derivative with respect to q^. 

8. Self-consistent equation for the Yukawa-type vertex T s (p,q). 

9. Another form for the second derivative of the auxiliary field propagator. 

10. Four-point fermion Green function: (a) a-exchange diagram, (b) pure ladder 
diagram. 

11. (ipijj) inserted Green function. 

12. RG functions in the symmetric and the M-dependent renormalizations: (a) (3 
function, (b) anomalous dimension. Solid line is for the symmetric renormaliza- 
tion (M = 0), while dotted line is for M = M d . 
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